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ABSTRACT 

Because  breaking  internal  waves  produces  most  of  the  turbulence  in  the  thermocline.  the  statistics  of  i.  the 
rate  of  turbulent  dissipation,  cannot  be  understood  apart  from  the  statistics  of  internal  wave  shear.  The  statistics 
of  e  and  shear  are  compared  for  two  sets  of  profiles  from  the  northeast  Pacific.  One  set.  PATCHEX,  has  internal 
wave  shear  close  to  the  Garrett  and  Munk  model,  but  the  other  set.  PATCHEX  north,  has  average  10-m  shear 
squared,  {S'?o),  about  four  times  larger  than  the  model. 

The  10-m  shear  components,  S^  and  S,,  were  measured  between  I  and  9  MPa  and  referenced  to  a  common 
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stratification  by  WKB  scaling.  The  scaled  components.  .S,  and  S,..  are  found  to  be  independent  and  normally 
distributed  with  zero  means,  as  assumed  by  Garrett  and  Munk.  This  readily  leads  to  analytic  forms  for  the 
probability  densities  of  S^g  and  Sto-  The  observed  probability  densities  of  Sjo  and  S*a  are  close  to  the  predicted 
forms,  and  both  are  strongly  skewed.  Moreover,  and  <riasto  constants,  independent  of  the  standard 
deviations  of  S,  and  5,..  The  probability  density  of  the  inverse  Richardson  number,  Ri,'o  *  S\o/(hi^).  is  a 
scaled  version  of  the  probability  density  of  55o-  The  PATCHEX  distnbuiion  cuts  off  near  Rim  =  4.  as  found 
by  Eriksen,  but  the  PATCHEX  north  distribution  extends  to  higher  values,  as  predicted  analytically.  Consequently  , 
a  cutoff  at  Ri  7o  =  4  is  not  a  universal  constraint. 

Over  depths  where  (A’^)  is  nearly  uniform,  the  probability  density  of0.5-m  <  can  be  approximated,  to  varying 
d^ees  of  accuracy,  as  the  sum  of  a  noise  variate  with  an  empirically  determined  distribution  and  a  lognormally 
distributed  variate  whose  parameters  can  be  estimated  using  a  minimum  chi-square  fitting  procedure.  The  0.5- 
m  t,  however,  are  far  from  being  uncorrelated,  a  circumstance  not  considered  by  Baker  and  Gibson  in  their 
analysis  of  microstructure  statistics.  Obtaining  approximately  uncorrelated  samples  requires  averaging  over  10 
m  for  PATCHEX  and  15  m  for  PATCHEX  north.  These  lengths  correspond  approximately  to  reciprocals  of 
the  wavenumbers  at  which  the  respective  shear  spectra  roll  off.  After  correcting  the  uncorrelated  t  samples  for 
noise  and  then  scaling  to  remove  the  dependence  on  stratification,  the  scaled  dissipation  rates. «  ■  <(  Ao/<  A'  > ), 
are  lognormally  distributed.  (Without  noise  correction  and  {N'‘)  scaling  the  dau  are  not  lognormal,  e.g.,  noise 
correction  and  scaling  with  {N')  and  <  A'*'^)  do  not  produce  lognoimality.) 

It  is  hypothesized  that  the  approximate  lognormality  of  bulk  ensembles  oft  results  from  generation  of  turbulence 
in  proportion  to  S*o-  Lognormality  is  well  established  for  isotropic  homogeneous  turbulence  (Gurvich  and 
Yaglom),  and  Yamazaki  and  Lueck  show  that  it  also  occurs  within  individual  turbulent  patches.  Bulk  ensembles 
from  the  thermocline,  however,  include  samples  from  many  sections  lacking  turbulence  as  well  as  from  a  wide 
range  of  uncorrelated  turbulent  events  at  different  evolutionary  stages.  Consequently,  the  bulk  data  do  not  meet 
the  criteria  used  to  demonstrate  lognormality  under  more  restricted  conditions.  If  the  authors  are  correct,  the 
high-amplitude  portion  of  <Af^)-scaled  bulk  ensembles  is  lognormal  or  nearly  so  owing  to  generation  of  the 
turbulence  by  a  highly  skewed  shear  moment.  As  another  consequence,  =  2.57  should  be  an  upper  bound 
for  10  m  ff|„;  when  the  turbulence  is  produced  by  the  breaking  of  random  internal  waves.  Because  man)(  ptarts 
of  the  profile  lack  turbulence,  sensor  noise  limits  the  «  distribution  to  smaller  spreads  than  those  of  Sio-  In 
practice  we  observe  ai„-,  =  1.2  when  Sjo  equals  GM76,  and  5^9  =  15  when  S|o  is  about  three  times  GM76. 
For  the  larger  spread,  95%  confidence  limits  requite  n  =»  60  for  accuracies  of  ±100%.  n  =  140  for  ±50%.  and 
n  =  20(X)  for  ±10%.  Owing  to  instrumental  uncertainties  in  <  estimates,  the  authors  suggest  accepting  less 
restrictive  confidence  limits  at  one  site  and  sampling  at  multiple  sites  to  estimate  average  dissipation  rates  in 
the  thermocline. 
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1.  Introduction 

For  the  past  20  years,  attempts  to  estimate  turbulence 
and  vertical  diffusivity  in  the  thermocline  have  been 
statistical  in  nature,  relying  on  bulk  averages  of  micro¬ 
structure  parameters  rather  than  on  detailed  time  series 
of  mixing  events.  At  the  same  time,  Garrett  and  Munk 
(1972a)  parameterized  the  average  characteristics  of 
the  background  internal  wave  field  in  terms  of  the  local 
buoyancy  frequency,  N.  Garrett  and  Munk  ( 1972b) 
then  applied  this  parameterization  to  the  rate  at  which 
internal  waves  randomly  break.  Munk  ( 198! )  and  De- 


1 


1778  JOURNAl  or  PHYSIC 

saubies  and  Smith  (1982)  refined  these  approaches, 
and  McComas  and  Miiller  (  1981  )  and  Henyey  et  ai. 

( 1986)  extended  the  analysis  by  showing  that  the  rate 
at  which  interactions  within  the  internal  wave  field 
move  energy  toward  small  scales  and  ultimate  dissi¬ 
pation  depends  on  A'  and  E.  the  dimensionless  energy 
density  of  the  internal  wave  field.  Ihat  is. 

X  (  A'*)/;',  where  (  )  denotes  an  arithmetic  mean. 

If  this  line  of  development  is  correct,  the  average  dis¬ 
sipation  rate,  (r),  depends  only  on  these  average  pa¬ 
rameters  and  is  thus  independent  of  the  physical 
mechanisms  producing  the  dissipation. 

Because  free-fall  profilers  do  not  adequately  measure 
E,  which  is  dominated  by  very  low  wavenumbers, 
Grcgg(  1989)  recasts  the  McComas  and  Miillerf  1981  ) 
and  Henyey  et  al.  ( 1986)  results  in  terms  of  the  variance 
of  internal  wave  shear, 

^  7  X  10  (Wkg  '].  (1) 

Ad  0(j,M 

where  A’o  =  0.0052  s  '.  Sn,  is  the  fourth  moment  of 
the  10-m  shear  produced  by  random  internal  waves, 
and  .So  1/  is  the  corresponding  moment  for  the  Garrett 
and  Munk  ( 1975)  model.  (This  model  is  known  as 
GM76  when  modifications  by  Cairns  and  Williams 
( 1976)  are  included.]  The  revised  scaling  was  applied 
to  six  sets  of  profiles  of  instantaneous  shear  and  dis¬ 
sipation  from  sites  where  internal  waves  provide  most 
of  the  shear.  For  these  data,  the  largest  ratio  of  <<>/ 
f<A'")  is  28.1  times  the  smallest  (c  is  the  kinematic 
viscosity  and  A'-)  is  a  measure  of  the  intensity 

of  turbulence  in  stratifiecl  fluids).  Applying  ( 1 )  reduces 
this  range  to  a  factor  of  .1.  Therefore,  in  view  of  the 
limited  variability  of  (A'*)  in  the  thcrmocline,  vari¬ 
ability  in  shear  appears  to  be  the  principal  factor  caus¬ 
ing  variations  in  /  «>.  Wavenumber  spectra  of  this  shear 
variability  are  given  by  Gregg  et  al.  (  1992 ). 

To  examine  whether  ( 1 )  applies  to  individual  lO-m 
estimates,  we  use  two  sets  of  profiles  taken  with  the 
multiscale  profiler  (MSP)  between  1  and  9  MPa  ( 100 
to  900  m)  during  the  Patches  Experiment  ( PATCHEX ) 
in  the  eastern  North  Pacific.  The  first  dataset  ( Gregg 
and  Sanford  1988)  contains  27  profiles  from  a  diffu¬ 
sively  stable  thermocline  in  which  the  variance  of 
finescale  shear  is  close  to  GM76.  Several  days  after  the 
main  PATCHEX  observations,  five  profiles,  referred 
to  as  PATCHEX  north,  were  taken  below  a  week-old 
coastal  jet  off  northern  California.  The  shear  variance 
of  the  second  dataset  is  about  four  times  GM76.  Of 
the  six  datasets  used  by  Gregg  (  1989)  to  compare  < 
and  shear,  only  the  two  sets  analy/.ed  here.  PATCfiFX 
and  PATCHEX  north,  contain  simultaneous  mea¬ 
surements  of  (  and  shear  witii  internal  waves  at  average 
or  above-average  intensity. 

fo  treat  the  observations  in  each  dataset  as  single 
ensembles,  the  A  dependence  of  the  profiles  is  removed 
with 
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Thus,  f  is  a  simple  scaling  oflhc  fourth  shear  moment; 
that  is,  combining  (  I  )  and  (  2  )  vickls 

<■„.  (),9.s9.St,.  (4) 

where  the  constant  has  units  of  meters  squared  per 
second,  and  A  scaling  of  shear  is  defined  in  the  next 
section.  In  Fig.  1.  is  plotted  on  the  v  axes,  and  <  is 
on  the  f  axes.  As  is  evident  in  the  top  panel,  lO-m  pairs 
are  weakly  correlated  and  have  trends  beginning  near 
login  Rii(l  =  U.  Ncvenheles.s.  the  large  scatter  dem¬ 
onstrates  that  equality  between  <  and  dcxrs  not  apply 
to  individual  estimates  but  to  ascrages.  Averaging  o\er 
101)  m  greatly  reduces  the  scatter,  and  using  SOI)  m 
produces  a  tight  distribution  about  the  dashed  line  rep¬ 
resenting  the  model. 

Because  (  I  )  applies  only  to  averages,  we  treat  the 
scaling  as  a  relation  between  two  random  processes. 
Specifically,  we 

•  test  how  well  our  observations  can  be  modeled  by 
random  processes, 

•  examine  to  what  degree  probability  densities  oft 
and  arc  lelated, 

•  form  confidence  limits  for  these  variables,  and 

•  examine  what  the  confidence  limits  imply  about 
sampling. 

Gargett  ( 1 990 )  suggests  that  scalings  other  than  (  1  ) 
provide  equally  good  fits  to  data.  This  may  be  true, 
but.  for  the  100-m  and  800-m  averages  plotted  in  Fig, 
I.  two-sample  Kolmogoroff-Smirnoff  tests  show  no 
significant  differences  between  the  e  and  distribu¬ 
tions.  Consequently,  we  cannot  reject  the  hypothesis 
that  in  each  case  e  and  arc  identical  populations. 
Therefore,  whether  or  not  (  1 )  is  a  unique  scaling,  it 
adequately  relates  averages  of  dissipation  and  shear. 
Determining  which  scaling  is  best  will  require  I )  ad¬ 
ditional  simultaneous  measurements  of  shear.  .V.  and 
<;  and  2)  careful  examination  of  the  statistics.  We  ini¬ 
tiate  the  latter  in  this  paper  and  focus  on  the  averaging 
required  to  collapse  widely  scattered  pairs  of  <  and 

to  the  close  equality  found  between  their  averages. 

In  section  2  the  normality  of  the  shear  components 
is  tested  and  probability  densities  arc  developed  for 
Riio',  and  ,S'|ii.  After  a  brief  review  of  the  lognor- 
maliiy  off  at  the  beginning  of  section  .^,  we  test  the 
lognormality  of  our  data  and  then  compare  the  prob¬ 
ability  densities  of?  and  In  section  4  we  present 
confidence  limits  for  the  shear  moments  and  for  log- 
niirmal  t  data,  and  we  also  discuss  implications  for 
sampling.  We  conclude  in  section  5  with  a  summary 
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PKI.  1 ,  5>cattcrplots  oflogm?  versus  l(igin<,»  for  averages  over  10  m.  1(K)  ni.  and  XOO  m.  PA  K'tUX 
data  arc  shown  by  dots  in  the  lop  panel  and  stars  in  the  lower  panels,  PA'ICMPX  north  data  arc 
shown  by  open  circles.  The  heavy  dashed  line  is  login*  logi„<,„- 


and  discussion.  Throughout  the  paper  we  use  non- 
parametric  statistical  tests  wherever  possible.  Often 
several  tests  are  applied  to  examine  different  aspects  of 
the  data.  To  maintain  continuity  in  the  discussion  while 
being  thorough  with  the  statistics,  many  details  of  the 
tests  arc  in  the  appendixes. 

2.  Statistics  of  shear  and  its  moments 

Our  goal  in  this  section  is  to  investigate  how  well 
instantaneous  shear  and  its  moments  can  he  repre¬ 


sented  as  random  processes,  as  assumed  h\  CjM76.  This 
involves  first  demonstrating  uncorrelatedncss.  in  the 
vertical  and  pairwise  between  components.  Ihen  we 
can  form  moments  of  the  data  and  test  whether  the 
distributions  fit  analytic  models.  Because  tests  applied 
to  individual  profiles  are  not  rigorous,  owing  to  small 
sample  sizes,  we  use  several  methods  for  each  hypoth¬ 
esis  we  test,  expecting  gross  departures  from  our  hy¬ 
potheses  to  be  apparent  in  all  tests.  We  note  minor 
depamires.  but  assume  that  their  presence  does  not 
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invalidate  further  development.  Details  of  statistical 
tests  used  in  this  section  are  in  appendices  A-C. 


a.  The  shear  components,  5,  and  Sy 

With  eastward  and  northward  flow  represented  by 
U  and  V,  the  shear  components  are  computed  as 

with  Az  =  10  m  (Gregg  and  Sanford  1988).  The  factor 
ofVz  .  1 1  corrects  for  the  attenuation  of  the  first-differ¬ 
ence  filter,  making  the  shear  variance  comparable  to 
an  integral  of  GM76  to  0.6  rad  m“‘  (0. 1  cycle  per  me¬ 
ter).  Systematic  changes  resulting  from  changes  in 
stratification  are  removed  with 


\ 

(Nyr., 


(6) 


Based  on  two  nonparametric  statistical  tests,  the  run 
test  and  the  spectral  distribution  function,  we  cannot 
reject  at  the  5%  level  of  significance  the  hypothesis  that 
for  PATCHEX  5*  and  Sy  are  vertically  independent, 
although  there  is  some  evidence  of  weak  vertical  cor¬ 
relation.  ( See  appendix  A  for  details.)  Neither  can  we 
reject  the  hypothesis  that  the  shear  components  are 
pairwise  independent  of  each  other.  We  thus  assume 
the  PAICHEX  shear  components  are  independent 
both  in  the  vertical  and  pairwise.  This,  however,  is  not 
true  for  PATCHEX  north,  as  one  of  the  five  profiles 
is  rejected  by  both  tests  (appendix  A). 

Means  and  standard  deviations  of  iV-scaled  and 
unsealed  shear  components  are  given  in  Table  1 .  Al¬ 
though  A  scaling  has  little  effect  on  the  standard  de¬ 
viations  for  PATCHEX,  it  increases  by  50%  the  stan¬ 
dard  deviations  for  PATCHEX  north.  Taking  cr^  and 
ffy  as  the  sample  standard  deviations  of  Sy  and  Sy,  we 
compute  normal  probability  densities  having  zero 
means  and  =  ( l/2)(ffy  +  ay).  (For  PATCHEX  5, 


Table  I.  Means  and  standard  deviations  of  10-m  shear  compo¬ 
nents.  The  PATCHEX  distributions  have  2187  samples,  and  the 
PATCHEX  north  distributions  have  405  samples.  Parameters  scaled 
with  are  overstruck  with  a  circumflex. 


Mean 

Standard  deviation 

PATCHEX 

Sy 

-2.20  X  10"’ 

2.50  X  10  ’ 

Sy 

-4.26  X  10"’ 

2.63  X  10  ’ 

Sy 

-4.88  X  10"’ 

2.74  X  10  ’ 

Sy 

-3.68  X  10"’ 

2.84  X  10"' 

PATCHEX  north 

Sy 

-2.68  X  10  ’ 

.3.33  X  10  ’ 

Sy 

-5.79  X  i0  ‘ 

3.20  X  10 

Sy 

-7.25  X  10  ’ 

4.91  X  10  ’ 

Sy 

7.82  X  lO  * 

4.74  X  10  ’ 

=  2.79X  10- 's',  and  for  PATCHEX  north  a,  =  4.83 
X  10"^  s'  ‘.)  As  seen  in  the  lower  panels  of  Fig.  2, 
overlaying  these  normal  probability  densities  and  his¬ 
tograms  of  the  observations  reveals  good  fits. 

Working  with  our  assumption  of  independence,  we 
apply  the  Lilliefors  test  statistic  (appendix  A)  to  test 
the  hypothesis  that  Sx  and  Sy  are  normally  distributed 
with  zero  means  and  common  variances.  We  find  that, 
at  the  5%  level  of  significance,  we  cannot  reject  this 
hypothesis.  A  visual  equivalent  of  this  test  is  displayed 
in  the  upper  panels  of  Fig.  2.  which  show  quantile- 
quantile  (g-g)  plots  (Chambers  et  al.  1983),  with  ob¬ 
served  distributions  along  the  vertical  axes  and  fitted 
normal  distributions  along  the  horizontal  axes.  If  the 
observed  distributions  were  identical  to  the  fitted  dis¬ 
tributions,  the  data  would  plot  along  diagonals  from 
lower  left  to  upper  right.  The  Lilliefors  test  statistics 
provide  a  way  of  plotting  95%  confidence  limits  to  as¬ 
sess  whether  ooserved  departures  from  the  diagonal 
are  consistent  with  random  deviations.  For  both  da¬ 
tasets,  Sy  remains  within  the  95%  confidence  limits,  as 
is  the  case  for  S,  ( not  shown ) .  Therefore,  we  can  treat 
the  shear  components  as  normally  distributed  with  zero 
means.  Under  this  assumption,  we  use  an  Ftest  to  find 
that,  at  the  5%  level  of  significance,  the  variances  of  Sy 
and  Sy  can  be  taken  to  be  the  same. 


b.  The  second  moment,  Sw.  and  the  inverse 
Richardson  number 


Unsealed  and  scaled  second  moments  are 


+  +  (7) 

When  Sy  and  Sy  are  normally  distributed  and  uncor¬ 
related,  the  probability  density  of  S'fo  can  be  calculated 
readily  because  it  is  the  sum  of  the  squares  of  Sy  and 
Sy.  The  result  is  an  exponential  distribution  (Papoulis 
1984), 


e\\>(-yl2a]) 

2o] 


H{y),  with  y 


Sw 


Si  +  S?; 
(8) 


//(>’)  =  1  for  ,v  >  0  and  //(,v)  =  0  for  v  <  0.  ( If  nor¬ 
malized  to  (t|  =  1 ,  Fy  is  X  ’  with  two  degrees  of  freedom.) 
The  g-g  plots  of  logioSto  show  the  PATCHEX  distri¬ 
bution  slightly  crossing  the  95%  confidence  limits  but 
all  the  PATCHEX  north  distribution  remaining  inside 
(Fig.  3).  (The  probability  density  of  logio^'Ki  is  derived 
in  appendix  B.) 

To  the  degree  that  the  gradient  Richardson  number 
is  not  influenced  by  fluctuations  in  induced  by  in¬ 
ternal  waves,  the  probability  density  of  logio  Ri  To  is 
simply  a  shifted  version  of  the  probability  density  of 
the  second  moment  of  the  scaled  shear;  that  is,  logm 
Rifo  =  logio^io  +  4.57.  where  logioAo  =  -4.57.  The 
critical  value  for  dynamic  instability  of  parallel  shear 
flows  is  Rir'  =4(Miles  196 1 ),  corresponding  to  logm 
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-0.02  0.00  0.02  -0.02  0.00  0.02 
S,  /  S’’  Sx  /  S’’ 

Fig.  2.  In  the  lower  panels,  normal  probability  densities  having  zero  means  and  standard  deviations  of 
a,  are  shown  by  heavy  lines.  The  observed  probability  densities  of  are  shaded.  In  the  upper  panels,  q~q 
plots  have  the  observed  (empirical)  distributions  along  the  vertical  axes  and  the  normal  (theoretical)  dis¬ 
tributions  along  horizontal  axes.  The  observed  distributions  fall  within  the  95%  confidence  limits  of  the 
theoretical  distributions  (dashed  lines),  which  are  much  wider  for  PATCHEX  north,  405  samples,  than  for 
PATCHEX,  2152  samples.  Appendix  A  describes  calculation  of  the  95%  confidence  limits. 


Rij'  =  0.6.  From  the  logio  RiTo  scale  in  Fig.  3,  we  see 
that  the  PATCHEX  distribution  drops  to  zero  just  to 
the  left  of  0.6,  similar  to  previous  reports  of  cutoffs  of 
Richardson  number  at  the  critical  value  ( Eriksen  1978, 
1982).  Only  0.2%  of  the  PATCHEX  samples  have  logm 
Ri7o  >  0.6,  compared  with  17.5%  exceeding  logm 
Ri  10  >  0.  The  observed  probability  density  for 
PATCHEX  north  also  follows  the  theoretical  predic¬ 
tion,  but  it  does  not  cut  off  at  the  critical  value;  9%  of 
the  samples  exceed  logic  RiFo  =  0.6,  and  56%  have 
logic  Ri7o  >0.  Replotting  in  more  conventional  form 
( Fig.  4)  shows  the  PATCHEX  north  Richardson  num¬ 
bers  peaking  sharply  between  0.25  and  0.5. 

Munk  (1981)  argues  that  saturation  of  the  shear 
spectrum  prevents  Richardson  numbers  from  dropping 
below  1  /4,  even  when  the  shear  rises  above  the  back¬ 
ground  state  modeled  in  GM76.  He  bases  his  argument 
on  observations  by  Eriksen  ( 1978 ),  who  computed  the 
Richardson  number  by  first-differencing  data  from 
current  meters  moored  7  m  apart  in  the  vertical.  Er¬ 
iksen  did  not  discuss  the  energy  level  of  the  internal 


waves,  but  review  of  his  data  shows  that  the  energy 
level  cannot  be  distinguished  from  GM76  (Eriksen 
1991,  personal  communication).  Thus,  Eriksen ’s  ob¬ 
servations  are  fully  consistent  with  our  distribution  for 
PATCHEX,  but  the  PATCHEX  north  distribution 
demonstrates  that  Eriksen’s  result  cannot  be  general¬ 
ized;  10-m  Richardson  numbers  smaller  than  one- 
fourth  occur  for  shears  above  GM76. 

c.  The  fourth  moment,  S/o 
Unsealed  and  A-scaled  fourth  moments  are 

5?o*(5?o)',  5to^(5fo)'.  (9) 

As  shown  in  appendix  B,  we  can  use  ( 8 )  to  obtain  the 
probability  density  of  namely, 

P  n(Fi)  exp(->>/2<r^) 

'  w  I  r/  \ I  xi  2  ^iy)* 

W  =  y'  =  Sto.  ( 10) 

Provided  that  and  Sy  are  normally  distributed  and 
uncorrelated,  the  standard  deviation  of  the  natural  log- 
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Fig.  3.  In  the  lower  panels,  observed  probability  densities  of  logio^io  (shaded)  are  compared  with  the 
predicted  form  (Bl ),  shown  as  heavy  curves.  The  light  curves  are  the  cumulative  distribution  functions  of 
the  observed  distributions.  Vertical  dotted  lines  mark  constant  values  of  logio  Ri  k'  .  The  critical  Richardson 
number  of  0.25  corresponds  to  log.o  Ri  =  0.6.  In  the  upper  panels,  q-q  plots  have  the  observed  ( empirical ) 
distributions  along  the  vertical  axes  and  the  theoretical  distributions  along  horizontal  axes.  Part  of  the 
PATCHEX  distribution  falls  slightly  outside  the  95%  confidence  limits,  but  none  of  the  PATCHEX  north 
distribution  exceeds  the  confidence  limits.  Confidence  limits  were  computed  using  the  exponential  case  of 
Mason  and  Bell  ( 1986). 


arithm  of  the  Sto,  finito’  ^  constant  equal  to  2.57 
and  does  not  change  with  variations  in  E.  The  distri¬ 
bution  of  logio^to  (Bl  1 )  is  highly  skewed,  with  a  long 
tail  extending  to  very  small  values.  Comparison  with 
the  observed  probability  densities  shows  parts  of  the 
PATCHEX  distribution  slightly  outside  the  95%  con¬ 
fidence  limits  on  q-q  plots  and  all  of  the  PATCHEX 
north  distribution  inside  the  limits  (Fig.  5). 

To  investigate  why  the  PATCHEX  distributions  fall 
just  outside  the  95%  confidence  limits,  we  examined 
the  vertical  correlation  structure  of  the  shear  moments. 
As  discussed  in  appendix  C,  we  find  no  evidence  of 
vertical  correlatedness  in  Slo,  but  there  are  indications 
of  temporal  correlatedness.  For  PATCHEX.  direct  ev¬ 
idence  comes  from  averaging  Sjo  and  5u)  vertically 
over  100  m  and  applying  run  tests  to  time  series  in 
each  pressure  range.  For  both  moments,  two  of  the 
eight  series  fail  the  run  test,  although  Sy  and  5,  pass. 
We  are  uncertain  of  the  reason,  but  believe  it  likely 
that  near-inertial  motions  dominate  the  shear.  Owing 


to  their  rotation,  they  rapidly  decorrelate  5,  and  S,  but 
not  Sw  and  5^)-  In  addition,  when  data  from  all  profiles 
in  the  same  cruise  are  treated  in  a  single  ensemble,  the 
standard  deviations  of  logio^fo  and  logio^to  do  not  de¬ 
crease  as  rapidly  as  they  should  for  independent  sam¬ 
ples  (appendix  C).  In  some  cases,  observed  standard 
deviations  exceed  predictions  by  more  than  50%.  The 
discrepancies  are  largest  when  averages  of  full  profiles 
are  compared,  which  suggests  weak  temporal  coher¬ 
ence. 

d  Summary 

Probability  densities  of  X,  and  X,  are  adequately  ap¬ 
proximated  with  normal  distributions  having  zero 
means  and  common  variance.  Empirical  tests  suggest 
that  S\  and  ,SV  are  uncorrelated  for  PATCHEX.  but 
PATCHEX  north  may  have  some  correlatedness.  .As¬ 
suming  uncorrclatcdness,  we  derive  analytic  forms  for 
the  probability  densities  of  logio-Vin  and  logioX^o.  Ob- 
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Fig  .  4.  Observed  probability  densities  of  Ri  m  ( shaded )  compared 
with  predicted  forms  (solid  lines).  At  low  values,  both  observed  dis- 
tnbutions  cut  off  as  predicted  theoretically.  For  P.ATCHE.X.  the  cutoff 
is  close  to  Riio  =  1/4.  For  PATCHEX  north,  the  cutoff  occurs  at 
smaller  values,  and  9%  of  its  samples  are  smaller  than  I  /4. 


served  probability  densities  of  these  moments  are  close 
to  the  predictions,  with  the  PATCHEX  north  distri¬ 
bution  within  the  95%  confidence  limits  and  the 
PATCHEX  distribution  slightly  outside  at  some  places. 
Because  this  pattern  is  contrary  to  what  would  be  ex¬ 
pected  from  the  correlation  tests  on  5,  and  S,  .  we  sus¬ 
pect  that  the  discrepancy  results  from  weak  temporal 
coherence  between  the  PATCHEX  profiles  and  that 
this  coherence  results  from  near-inertial  motions. 
Probability  densities  of  Richardson  number,  computed 
using  the  mean  stratification,  agree  with  the  predicted 
distributions.  The  PATCHEX  distribution  cuts  olf  near 
Riio  =  1/4.  as  found  by  Eriksen  (1978),  but  the 
PATCHEX  north  distribution  extends  to  smaller  val¬ 
ues,  Therefore,  a  cutoff  at  Ri  n,  =  1  /4  is  not  a  universal 
condition, 

3.  Statistics  of  dissipation 

Turbulence  is  inherently  intermittent.  Even  when 
turbulence  is  statistically  homogeneous,  local  concen¬ 
trations  of  small-scale  shear  determine  average  dissi¬ 
pation  rates,  (<),  throughout  the  entire  volume.  Con¬ 


sequently,  the  sampling  problem  is  to  observe  enough 
of  the  infrequent  large  values  to  ensure  precise  averages. 

Tor  statistically  homogeneous  and  isotropic  turbu¬ 
lent  flows.  Gurvich  and  Yaglom  (  1967  )  show  that  the 
turbulent  cascade  leads  to  lognormal  probability  den¬ 
sities  of  t:  that  is.  Ine  is  normally  distributed.  Lognormal 
distributions  are  typical  for  positive-definite  variables 
produced  by  repeated  multiplicative  applications  of  a 
random  process.  For  example,  the  sizes  of  sand  grains 
broken  off  nearby  rocks  are  lognormal,  as  are  dissi¬ 
pation  rates  of  the  smallest  eddies  in  a  turbulent  cas¬ 
cade.  Lognormal  probability  densities  are  parameter¬ 
ized  by  /xin,  and  crin, .  the  mean  and  standard  deviation 
of  ln«.  Gurvich  and  Ya^om  ( 1967)  argue  that  o,n. . 
which  measures  intermittencc.  and  hence  the  difficulty 
of  sampling,  increases  with  Reynolds  number.  At  pres¬ 
ent,  however,  there  is  no  way  to  estimate  crin,  a  priori, 
even  in  homogeneous  turbulence  (Yamazaki  and 
Lueck  1990). 

The  airfoil  probes  used  on  MSP  to  detect  velocity 
microstructure  also  respond  to  thermal  transients  hav¬ 
ing  frequencies  of  0.1-1  Hz  (Osborn  and  Crawford 
1980).  At  the  MSP  fall  rate,  these  frequencies  corre¬ 
spond  to  wavenumbers  of  0.3-3  cpm.  To  minimize 
the  effect  of  thermal  transients  and  to  detect  short- 
length  changes  in  turbulent  intensity,  we  compute  e  by 
integrating  spectra  taken  over  0.5-m  data  windows. 

In  examining  the  statistics  of  e.  we  must  consider 
w'hether  the  samples  are  correlated.  Only  with  uncor¬ 
related  samples  can  we  accurately  bound  the  degree  of 
intermittencc  and  estimate  confidence  limits  for  bulk 
ensembles  of  dissipation.  Uncorrelatedness  has  not 
been  addressed  previously — investigators  have  merely 
assumed  it.  regardless  of  sample  length  or  local  cir¬ 
cumstances.  For  example.  Baker  and  Gibson  ( 1987) 
assume  uncorrelatedness  for  0.23-1.15-m  samples 
along  horizontal  tows  taken  by  Washburn  and  Gibson 
( 1984)  and  for  100-200-m  samples  along  profiles  taken 
by  Gregg  ( 1977).  Confidence  limits  computed  using 
standard  deviations  from  positively  correlated  data  are 
usually  optimistic:  that  is.  they  give  tighter  limits  than 
warranted. 

a.  Prohahilily  densities  oj  0.5-m  t 

Probability  histograms  of  logiof  for  the  0.5-m  sam¬ 
ples  are  highly  skewed,  with  long  tails  al  large  magni¬ 
tudes  and  sharp  cutoffs  al  low  magnitudes  (Fig.  6). 
Probability  distributions  have  long  been  described  as 
resulting  from  approximately  lognormal  oceanic  signals 
at  high  magnitudes  and  noise  at  low  magnitudes.  Ob¬ 
servations  are  considered  noise  when  their  spectral  level 
drops  to  a  low  level  that  is  only  approximately  constant 
because  the  ensembles  are  collected  using  probes  having 
difl'erent  noise  levels  and  sensitivities.  For  MSP  this 
level  is  10  W  kg  ‘  (Gregg  and  Sanford  1988).  Be¬ 
cause  (  is  positive  definite,  noise  distributions  arc  not 
expected  to  be  normal  but  are  closer  to  lognormal. 
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Fig.  S.  In  the  lower  panels,  the  observed  distributions  oflogio.$io  (shaded)  agree  well  with  the  prediction 
(B1 1 ).  In  the  upper  panels,  q-q  plo'i  based  on  modified  Lilliefors  statistics  show  parts  of  the  PATCHEX 
distribution  slightly  exceeding  the  95%  confidence  limits.  None  of  the  PATCHEX  north  distribution  falls 
beyond  the  bounds. 


Distributions  that  are  mostly  noise  were  observed  by 
MSP  on  the  equator  during  TROPIC  HEAT  2  in  April 
1987.  Below  the  undercurrent,  zones  several  hundred 
meters  thick  frequently  contain  little  or  no  measurable 
dissipation.  In  Fig.  6  the  distribution  of «  from  9  to  9.5 
MPa  during  TROPIC  HEAT  2  is  seen  to  be  similar  to 
the  low-magnitude  portion  of  the  PATCHEX  distri¬ 
bution. 

We  can  estimate  the  lognormal  parameters  by  fitting 
the  observed  histograms  to  a  function  that  is  the  sum 
of  the  TROPIC  HEAT  2  noise  distribution  and  a  log¬ 
normal  distribution.  Taking  data  in  increments  of  1 
MPa,  we  use  a  minimum  procedure  to  estimate 
means,  ,  and  standard  deviations,  <T|n< .  These  esti¬ 
mates  are  given  as  ^i„,  and  Sin,.  As  described  in  appen¬ 
dix  D,  the  procedure  works  best  for  PATCHEX  north, 
which  is  less  affected  by  noise.  For  PATCHEX, 
decreases  strongly  with  pressure,  and  increases 
gradually  from  ffi„,  =  2.9  at  2-3  MPa  to  =  3.2  at 
8-9  MPa  (Fig.  7).  Its  average  is  3.01  ±  0.13.  For 
PATCHEX  north,  jl\„,  increases  only  between  2  and  5 
MPa,  reflecting  the  nonuniform  vertical  distribution 


of  dissipation  and  shear,  and  ciin,  shows  no  trend  and 
has  an  average  of  2.56  ±0.17. 

b.  Averaging  0.5-m  t  lo  obtain  uncorrelated 

estimates 

Even  casual  inspection  of  profiles  of  0.5-m  (  shows 
vertical  trends,  with  patches  of  contiguous  large  values 
separated  by  zones  at  the  noise  level  (Gregget  al.  1986). 
Consequently,  applying  the  run  test  to  unsealed  e  be¬ 
tween  2.5  and  5.0  MPa,  where  N  is  uniform,  results  in 
all  profiles  being  rejected  (Table  2).  Increasing  the  av¬ 
eraging  length,  Z.;,  slowly  decreases  rejections.  For 
PATCHEX,  T;  =  8  m  yields  3  rejections  in  27  profiles, 
a  ratio  found  previously  for  random  rejections  of  un¬ 
correlated  data  (appendix  A).  Further  increasing  to 
10  m  removes  all  rejections,  and  we  treat  !0-m 
PATCHEX  averages  as  uncorrelated.  PATCHEX 
north  requires  averaging  over  10  lo  15m  for  no  rejec¬ 
tions.  In  both  cases,  the  minimum  averaging  length 
required  for  uncorrelated  samples  corresponds  to  the 
reciprocal  of  the  wavenumber  (in  cycles  per  meter)  at 
which  the  shear  spectra  roll  off  (Gregg  et  al.  1992). 
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Fig  .  6.  Probability  densities  of  all  0.5-tn  logio<  samples  between  I 
and  9  MPa  for  PATCHEX  (shaded)  and  for  PATCHEX  north  (solid 
line).  They  are  compared  with  a  probability  density  obtained  between 
9  and  9.5  MPa  on  the  equator  during  TROPIC  HEAT  2  (dashed 
line).  The  latter  appears  to  be  nearly  all  noise. 


c.  Noise  correction  of  the  IQ-m  estimates 

To  facilitate  comparisons  of  the  datasets,  we  average 
€  over  10  m  and  scale  with  to  form  e  ( 3).  Testing 
these  profiles  for  uncorrelatedness  over  the  full  range 
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of  1-9  MPa,  1 1  PATCHEX  records  fail  the  run  test 
and  6  fail  with  the  spectral  distribution  function  (Table 
3).  These  failures  result  from  the  increasing  etfect  of 
noise  below  5  MPa.  To  correct  for  the  noise,  we  set  to 
zero  all  0.5-m  samples  with  t  <  10  Out  of  2187 
lO-m  intervals  for  PATCHEX,  35  have  no  0.5-m  sam¬ 
ples  with  f  >  10  Omitting  these,  only  1  of  the 
PATCHEX  profiles  fails  the  run  test,  but  6  still  tail 
with  the  spectral  distribution  function  (Table  3). 
Therefore,  a  minority  of  the  records  have  some  vcnical 
correlatedness,  even  after  noise  correction.  Because 
noise  is  less  important  for  PATCHEX  north,  no  lO-m 
intervals  had  to  be  dropped,  and  noise  correction  does 
not  affect  the  run  test  or  the  spectral  distribution  func¬ 
tion  (Table  3).  Removing  noisy  values  decreases  the 
magnitudes  of  the  smallest  10-m  averages,  thereby  in¬ 
creasing  ffin,  for  both  datasets  (Table  4). 

d.  Probability  densities  of  noise-corrected  10-m  t 

Applying  (N^)  scaling  to  the  noise-corrected  data 
affects  the  standard  deviation  differently  for  the  two 
records  because  they  have  differing  vertical  distri¬ 
butions  of  dissipation  (Table  4).  Taking  the  (N^)- 
scaled  and  noise-corrected  distributions  as  the  most 
accurate,  for  PATCHEX  the  sample  mean  of  Int, 
Min!  -  -22.240  and  the  sample  standard  deviation  of 
Inl,  a,„;  =  1.188;  for  PATCHEX  north  Mm;  =  -20.270 
and  =  1 .463. 

After  noise  correction  and  (N^)  scaling,  both  da¬ 
tasets  pass  Q-q  tests  for  lognormality  (Fig.  8).  With 
only  scaling  or  with  only  noise  correction,  neither 
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Fig  .  7.  Parameters  for  0.5-m  ln«,  obtained  by  X  ^  fits  to  data  taken 
in  1  MPa  increments.  PATCHEX  values  are  shown  by  open  circles 
and  PATCHEX  north  by  crosses.  Reflecting  different  vertical  distri¬ 
butions  of  (f),  iiia.  decreases  rapidly  with  pressure  for  PATCHEX 
but  increases  slightly  for  PATCHEX  north.  For  both  data-sets,  esti¬ 
mates  of  the  standard  deviation,  are  nearly  constant  below  the 
near-surface  high- Af  duct. 


Table  2.  Run  tests  on  averages  of  0.5-m  <  between  2.5  and  5.0 
MPa.  Because  noise  contamination  is  insignificant  and  the  stratili- 
cation  nearly  uniform,  no  corrections  arc  applied  to  «.  The  number 
f,f  rr;-r«'ons  is  left  o*"  ’he  solidus,  and  the  number  of  samples  is  right. 
The  number  of  samples  in  each  prolile  decreases  with  increasing 
averaging  length.  .  All  profiles  fail  for  =  0.5  m,  which  represents 
no  averaging;  for  example,  for  PATCHEX  there  are  27  rejections  for 
27  profiles.  Rejections  decrease  with  averaging;  there  are  none  for 
=  10  m  for  PATCHEX  and  -  15  m  for  PATCHEX  north. 
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I  he  number  of  rejections  is  to  the  left  and  the  numtx-r  of  profiles  is 
to  the  right. 
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distribution  passes,  owing  to  vertical  trends  remaining 
in  the  profiles.  We  also  scaled  with  A  '  and  A  '  pro¬ 
posed  by  Gargett  (1990)  to  remove  vertical  trends. 
These  scalings  do  not  completely  remove  vertical 
trends,  and  owing  to  the  trends  the  scaled  data  do  not 
pass  for  lognormal. 

Averaging  the  lO-m  samples  over  100  and  800  m 
reduces  ^in,  by  factors  of  2-.').  considerably  less  than 
the  reduction  expected  for  uncorrelated  samples  ( Table 
5).  For  example,  averaging  ten  of  the  PATCHEX  sam¬ 
ples  reduces  (T|„;  from  1.19  to  0.64  instead  of  to  0.47 
as  expected.  This  is  similar  lo  the  siower-than-cxpected 
decrease  in  ffin.vf,,. 

e.  Comparison  with  prohahiliiy  densities  of  S^j,, 

Scaling  the  5m  distributions  as  and  overlaying 
them  on  the  noise-corrected  <  distributions  reveal  that 
both  distributions  are  similar  at  large  magnitudes  ( Fig. 
9).  (Appendix  E  includes  a  more  quantitative  com¬ 
parison  using  ^-f/ plots.)  Agreement  at  high  magnitudes 
is  important  because  the  high-magnitude  tails  heavily 
influence  the  mean  values,  or  first  moments.  The  ex¬ 
pected  value  of  is  obtained  by  using  (4  )  and  (  B.1): 

=  0.959E[St„]  =  0,9.‘>9  X  (II) 


It  is  estimated  by  0.9si)  •  ?<« '.  I  he  expected  value  of 
(  is 

/■.  f  f  1  e\pl  p;,,,  ■  n.Sfi  ),  (12) 

It  IS  estimated  In  espi  ij„..  ■  i,,.  I.  i  oi  P  \  I  (  I II  A. 

the  estimaie  of /,  I  ('  )is4.6s  .  ip  '  aiul  the  esiimate 
ol  /  |i|  is  4.4.''  •  10  I  hese  are  witliin  2  ot 
4.55  '  10  ".  the  average  of  ihe  noise-coi  rected  aiut 
A'  -sealed  dissi|xiiion  data.  1  or  P\  l<  111  .\  north, 
the  estimate  ot  /.  [  i  ,  |  is  4. !  S  •  ]0  '  and  the  estimate 
of  /.[()  is  4  59  ■  10  which  are  within  20'-  ol  i 
5,1,5  -  10 

Because  i  is  much  more  all'ectcd  by  noise  than  is 
f  ,  their  probability  densities  do  not  agree  at  low  mag¬ 
nitudes,  particularly  for  PATCTIEX  ( Fig.  9 ).  1  o  dem 
onstrale  the  elfect  of  noise,  vve  crudely  simulated  it  by 
adding  5  X  10  "(.\(,/.V)-  to  the  10-m  5m  values  and 
then  converted  them  to  .  As  seen  in  the  lower  panels 
of  Fig.  9,  adding  the  noise  greatly  improves  agreement 
at  low  magnitudes. 

4.  GunHdcnce  limits  and  sampling 

a  I'he  seamd  oioment  ol  shear.  .Sj„ 

Because  .S'fo  is  a  rescaled  x'  random  variable  with 
two  degrees  of  freedom,  confidence  limits  for  /:'(.Vi,i) 
=  2(t:  can  readily  be  found  based  on  an  average  of  n 
independent  samples.  (.Shi;  =(l/n)I',  'i  .S' i,,.,.  These 
limits  can  be  stated  as  the  interval 

\L„  X  <  SjiiS.  (.  I,  X  (Sjo/].  (1.5) 

As  shown  in  appendix  F, 

,  ,  2n 

..s  unci  L  fj  2  *  (  1  ) 

'  ‘  2nJ  \  ‘  /'\/2  ' »  I 

where  .V is  the  r/th  percentage  point  of  the  X  ’  dis¬ 
tribution  with  V  degrees  of  freedom,  and  p  is  the  con¬ 
fidence  level.  For/7  ~  0.95,  averaging  10  samples  gives 
/.i,i  =  0.585  and  Gk,  =  2.09  (Fig.  10).  As  n  increases. 


Tabi  r  4.  Statistics  of  lO-m  dissipation  between  1  and  9  MPa.  f  or  both  datasets  the  lO-m  estimates  are  formed  in  three  wavs.  Anthmelic 
averages  include  all  20  of  the  ().5-m  samples  in  the  lO-m  inurval.  Noise-iorreeted  lO-m  values  are  eonaputed  by  setting  to  zero  0  5-m 
samples  less  than  10  W  kg  ‘.  Tor  PATClItX  of  the  lO-m  averages  contained  no  values  larger  than  the  noise  and  were  dropped  from 
the  cc'rrected  ensemble.  Tinally.  the  notse-correeted  estimates  are  scaled  using  (,f|. 


Data 

n 

t 

Min. 

PAK  III  X 

arithmetic 

2IS7 

4. 11 

■  10 

22, .^5' 

I.Ml 

noise  corrected 

2152 

4.02 

10 

22.700 

I..V-A 

noise  corrected  and  -  A  ’,,  sealed 

2I.X2 

4.5.'i 

■  10 

22.240 
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PA  TCHI  X  north 

simple  average 

405 

2.12 

'  III  ■' 

20,d7X 

1.101 

noise  corrected 

405 

2  10 

r  10  ■* 

21.057 

i.4nx 

noise  corrected  and  A\  scaled 

405 

5. 1 1 

•  10  ^ 

2(f270 

1.401 

AK.ISI  IW 


( .  K  I  <  I  ( ;  I  I  \  1 


log,o/- 

l  ui.  ><.  Prohabilil>  densities  of  .V-  -sealed  and  noisc-correcicd  Hl-m  dissipation  rates  In  tne  lower 
panels,  observed  probabilits  densities  ( shaded  I  of  logi,i(  arc  compared  with  normal  distributions  basing  the 
same  means  and  standard  deviations  In  the  upper  panels.  <./-'/  plots  show  that  both  diaiibutions  fall  withm 
the  y.'s’-  confidence  limits.  Ib  construct  the  (/-</  plots,  we  use  the  mcxiilied  l.illiefors  test  with  case  ,t. 
estimating  means  and  standard  deviations  from  the  distribution. 


and  {  „  slowly  converge  lo  1;  averaging  the  27 
PATCHliX  profiles  between  500  and  600  m  gives  /.i^o 
=  0.89 1  and  t  =1.1 3.  and  averaging  all  PATC'HHX 
data  having  matching  noise  corrected  <  vields 
0.959  and  L  -,.:  1.04. 

h  Thi’  fourth  montent  of  shear.  .Sj/j 

Confidence  limits  for  /U.S’u,}  =  8(t?  2[/-.'(.Vi„))’ 

can  also  be  based  on  /  resulting  in 

[I.;,  X  2(<'.Vt„>)x  X  2(/.Vf,. .>)-].  (15) 

where  2((.S’T(,)  )’  is  an  estimate  of  (.Vu,  >  for  normally 
distributed  shear  components.  Alternatively,  confi¬ 
dence  lirnits  can  also  be  based  on  CS  ni  )  =  (1/n) 
X  II;  '|  .Vto,  and  take  the  form 

( X  ( .S  III  / .  {  ,1  X  f  .S  >  1  ■  (16) 

As  discussed  in  appendi.x  F.  /.;,  and  C;,  arcconsidcrabh 
more  difficult  to  obtain  than  /.„and  L  '„.  Moreover,  for 
n  >  I  the  expected  distance  between  the  upper  and 
lower  confidence  limits  based  on  (.S'm/  is  .3''<-12'f 


larger  than  the  expected  distance  based  on  f  ,S'],i  -.  For 
these  reasons,  the  confidence  limits  for  /•,(  .S't,,)  in  (  1 5  ) 
arc  preferable. 

Because  and  Cr,  depend  only  on  u.  they  are  the 
stimc  for  both  datasets.  Forming  the  l()()-m  averages 
plotted  in  F'ig.  I  requires  u  -  10,  resulting  in  /,  t,, 

().58s’  0.342  and  Cj,,  -  2.09'  4.36,  llsing(4', 

to  convert  .St,  to  f,« .  we  also  multiply  the  multiplicative 
confidence  limits  by  0.959.  Thus,  on  the  logarithmic 
axis  for  i,„  in  Fig.  1,  the  95';  confidence  limits  for  the 
100-m  values  correspond  to  0,48  and  i-0.62.  f  or  the 
80()-m  averages  in  F'ig.  1 .  the  limits  are  0.632  and  1 .52. 
corresponding  to  -  0.20  and  +0, 18  on  the  logarithmic 
axis. 

I'o  form  confidence  limits  for  the  expected  value, 
of  all  data  for  each  set,  for  PA  FCMFX  /.bs. 

0.920  and  (  tv:  1-08.  Because  the  estimate  of 

/.  I  «,„  ]  is  4.65  >  10  "’.the  95',  confidence  limits  are 
(4.IO  .V10  "',4.82  ■  10  "’).,Similarly.foi  [’AK'HFX 
north  (-io^  0.909  and  (  t,.  1 . 1  I .  With  an  estimate 

of  the  expected  value  equal  to  4,18  >  lo  gs'  ,  con¬ 
fidence  limits  are  (  3.64  x  10  4.45  >  10 
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TaHM:  Changes  in  jji„-  and  with  averaging  length.  F.sliinales  of  sample  mean'  and  slandard  dcMations  ol  In  aie  gnen  I",  f,  . 

and  a,„-.  lespceliveU.  .Averaging  lO-m  noise  eorreeled  i  over  ilK)  and  XtKl  m  invreases  and  deeieases  a;„,  I  or  vomiiansiin  ihe  'allies  m 
parentheses  are  the  ones  expected  for  averaging  lulls  indeix'iidenl  and  lognornialK  dislnhuted  m  in  vim()les  I  lies  ssere  lali  iil.ited  h.  .1 
Taylor  series  expansion  and  bs  Monte  Carlo  simulation 


PA  ICHFX  PA  It  HF\  north 
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f  The  dissipation  rate,  i 

For  lognormal  distributions,  the  maximum 
likelihood  estimate  of  /:'[?]  is  calculated  by  evaluating 
(12)  with  the  sample  mean  ^i,,;  and  sample  standard 
deviation  ain.  of  InF  In  appendix  G  we  discuss  some 
difficulties  with  Baker  and  Gibson's  (  1 987 )  use  of  this 
estimate. 

I-and  ( 1975 )  gives  tables  for  computing  exact  con¬ 
fidence  limits  for  /;'(?]  based  upon  mi„;  and  Sin,  .  These 
limits  can  be  written  as  (fmicT,,  emicTji]  where  I.„  and 


(  „  can  iHf  obtained  from  h:s  table.s.  Land  ;  1^72  )  dis- 
cu.sse.s  an  appruximation  to  these  limits  due  to  F).  R, 
Cox.  I \)r  95*', ■  confidence  Iimit.s  this  apprsiximalion 
takes  the  simple  form 

1  <nil.-  CXp(  -■  1  .96-)  ),  .'„„n  o.vpl  !  .90-y  )  j  .  (17) 

where 

7  ""  [  e'ncJn  f  nt„j2(n  4  |  (L  '  (  |  ,s  ) 

I  Baker  and  Gibson  (  1987)  obtained  a  similar  result). 
While  this  approximation  is  useful  for  large  sample 


■11  10  9  8  10  9  8 


log,o  e  log,(,  r 

F'lo.  9  Comparisons  of  probabilily  densities  of  l(l-m  logmi  (shaded)  and  lognu',,  (heavy  line)  The  i 
distributions  have  been  corrected  for  noise  and  scalerl  with  '  .V*  /.In  the  upper  p.inels  the  >  and  vlisinbuiions 
agree  much  better  at  high  magnitudes  than  at  lovx  magnitudes.  The  lovxcr  panels  repeat  the  comparison  after 
noise  similar  to  that  modeled  for  <  has  been  added  to  the  .  This  greatly  improves  the  overall  agreement 
between  the  <  and  f,«  probabilily  densities,  particularly  for  PA  fCIIFX 
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Fig.  10.  Multipliers  for  upper  and  lower  95%  confidence  limits  of 
averages  of  n  independent  samples  of  the  second  and  fourth  moments 
of  shear.  These  multipliers  are  t/„  and  L,  in  ( 1 3 )  for  5,0  and  Ul  and 
Li  in  (15)  for  5*0- 


sizes  and  moderate  values  of  cri„;,  it  can  give  large  errors 
for  small  sample  sizes.  For  example,  we  take  «  =  10 
with  cr,„;  =  1.19  for  PATCHEX  (Table  5).  Tht  ap¬ 
proximations  give  (0.39<mie,  2.57<mie)  compared  to 
(0.45?mie,  6.46emie)  for  the  exact  95%  confidence  limits 
[these  limits  correspond  to  (-0.35  4-  logiofmie^  +  0.81 
+  logio«mie)  on  logarithmic  axes]. 

The  multipliers  L„  and  U„  are  plotted  in  Fig.  1 1  for 
(Tin;  -  0.5,  1,  1.5,  2,  2.5,  and  3,  covering  the  range  of 
a\„i  likely  to  be  obtained  from  observations.  As  an  ex¬ 
ample,  the  values  in  the  first  row  of  Table  5  indicate 
that  =  exp(-22.24  +  0.5  X  1.19^)  =  4.45  X  10"'“ 
for  PATCHEX  and  «mie  =  exp(  -20.27  +  0.5  X  1.46^) 
=  4.57  X  10  ’  for  PATCHEX  north.  The  correspond¬ 
ing  95%  confidence  limits  for  £,'[«]  are  (4.17  X  10 
4.76  X  10"'“)  for  PATCHEX  (based  upon  «  =  2152 
samples)  and  (3.76  X  10  5.69  X  10  “)  for 

PATCHEX  north  (n  =  405  samples).  The  corre¬ 
sponding  multipliers  are  L2i^2  =  0.937  and  C2i,w  =  1  07 
for  PATCHEX  and  U),  =  0.823  and  =1.25  for 
PATCHEX  north.  For  both  datasets  the  confidence 
limits  include  the  estimates  obtained  using  irt  in  ( 1 1 ) 
and  the  computed  sample  averages  of  i. 

5.  Summary  and  discussion 

a.  Summary 

Statistics  of  shear  and  dissipation  are  compared  for 
two  sets  of  profiles  taken  when  the  internal  wave  field 


appeared  statistically  homogeneous  over  distances  of 
kilometers  and  steady  over  intervals  of  days.  For  one 
set,  PATCHEX.  the  internal  waves  are  close  to  GM76; 
for  the  other  set,  PATCHEX  north,  the  root-mean- 
square  10-m  shear  is  three  times  GM76.  Our  findings 
are  as  follows. 

1 )  After  scaling  the  shear  components.  5,  and 
with  ).  the  scaled  components,  and  5, ,  are  nor¬ 
mally  distributed  with  zero  means  and  are  uncorrelated 
with  each  other.  This  readily  leads  to  analytic  foirns 
for  the  probability  density  functions  of  Sio  and  Sto. 
The  latter  is  highly  skewed,  and  the  standard  deviation 
of  ln.9to  is  a  constant,  ^  2.57.  For  PATCHEX 
north  the  measured  distributions  of  Sw  and  S^o  are 
within  the  95%  confidence  limits  of  the  analytic  forms. 
For  PATCHEX  the  measured  distributions  are  similar 
to  the  analytical  forms,  but  they  slightly  exceed  the 
95%  limits. 

2)  For  both  datasets  the  probability  densities  of  the 

inverse  Richardson  number,  Rifo  =  are 

close  to  those  predicted.  Consistent  with  Eriksen 
(1978),  the  PATCHEX  distribution  cuts  off  near  4. 
but  9%  of  the  PATCHEX  north  distribution  has 
Ri7o  >  4.  Therefore,  the  cutoff  at  Ri^J  =  4  is  not  a 
universal  constraint. 

3)  For  bulk  ensembles  of  0.5-m  t  collected  where 
(N^)  is  nearly  constant,  the  lognormal  parameters. 
fi,„,  and  can  be  estimated  objectively  by  minimiz- 


Fio.  1 1.  Multipliers  for  upper  ( (  „}an<j  lower)  L„)9S't;  confidence 
limits  of  averages  of  lognormal  distributions,  as  functions  of  the 
number  of  samples,  n.  and  of  the  standard  deviation  of  the  natural 
logarithm  of  the  variable.  The  multipliers  arc  obtained  by  interpolating 
the  exact  confidence  limits  given  by  land  ( 1 975 )  and  are  giicn  for 
(Tin;  spanning  the  range  expected  in  the  ocean 
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ing  X"  fits  to  the  sum  of  a  lognormal  distribution  and 
an  empirically  determined  noise  distribution.  In  our 
cases,  the  noise  distribution  was  observed  in  an  un¬ 
usually  quiet  section  below  the  equatorial  undercurrent. 
Fits  to  0.5-m  e  taken  in  bins  of  1  MPa  vary  in  quality 
and  have  averages  of  irin,  equaling  .7.01  ±  0.1.7  for 
PATCHEX  and  2.56  ±  0. 17  for  PATCH  EX  north. 

4)  Half-meter  e  are  not  uncorrelated  but  must  be 
averaged  vertically  to  form  uncorrclated  samples,  fhe 
minimum  averaging  length  is  10  m  for  PATCHEX 
and  about  15  m  for  PATCHEX  north.  Both  averaging 
lengths  correspond  approximately  to  the  reciprocal  of 
the  wavenumber  at  which  their  respective  shear  spectra 
roll  off  ( Gregg  et  al.  1992). 

5)  After  correcting  for  noise  and  scaling  t  with 

(A"),  probability  densities  of  10-m  .V-scaled  dissipa¬ 
tion  rates,  t.  are  lognormal.  (Scaling  with  A'"'  and 
.'  A  '  -)  does  not  yield  lognormalit\.|  Thus,  the  loga¬ 
rithms  of  10-m  scaled  dissipation  arc  normal.  As  they 
are  also  uncorrelated,  they  can  be  assumed  to  be  in¬ 
dependent,  These  distributions  have  fT|„,-  1.19  for 

PATCHEX  and  1.46  for  PATCHEX  north,  fhe  10- 
m  samples,  however,  arc  not  fulls  independent; 
namely,  further  averaging  does  not  reduce  i7|,„  as  rap¬ 
idly  as  expected  for  independent  samples.  This  appears 
to  be  the  result  of  weak  vertical  correlation  within  pro¬ 
files  and  temporal  correlation  between  profiles. 

6)  For  10-m  samples,  the  probability  densities  of 
10-m  logt  and  log.S'fd  ( represented  as  loge,., )  arc  similar 
at  high  magnitudes,  and  mean  values  oft  and  differ 
by  only  2Fr-20'7 .  As  a  consequence  of  the  frequent 
lack  of  turbulence  at  low  shear  magnitudes,  noise  affects 
much  of  the  low-magnitude  probability  density  oft, 
causing  it  to  differ  from  the  t  distribution. 

7 )  Our  assumption  that  the  data  result  from  random 
processes  appears  valid  since  the  confidence  limits  cap¬ 
ture  the  collapse  with  averaging  ( Fig.  1 ) .  T  his  also  lends 
credence  to  our  approach  in  testing  for  independence. 

8 )  Ten  independent  measurements  arc  needed  to 
obtain  95'iV  confidence  limits  of  a  factor  of  2  for 
T(Ai(i)  =  2a;.  Owing  to  the  higher  skewness  of  its 
distribution,  similar  confidence  limits  for  .S'to  require 
as  few  as  40  independent  samples.  Sixty  independent 
dissipation  samples  are  also  needed  for  the  same  degree 
of  confidence  in  1 0-m  dissipation  rates  ( these  arc  com¬ 
puted  using  ain,  s;  1 .5  in  maximum  likelihood  estimates 
for  /•-'[?])•  fightening  the  confidence  limits  rapidly  in¬ 
creases  the  number  of  samples  required. 

h.  Discussion 

Wc  average  and  then  decimate  to  obtain  indepen¬ 
dent  samples.  As  a  reviewer  correctly  notes,  indepen¬ 
dent  samples  can  also  be  obtained  by  decimation  with¬ 
out  averaging.  Although  dccimation-t)nly  is  suitable 
for  some  studies,  we  need  to  average  and  decimate  to 
retain  the  average  dissipation  accompanying  the  inter¬ 


nal  wave  shear  over  the  same  vertical  interval.  Also, 
10-m  averaged  <  data  approximately  tolluw  a  lognormal 
distributimi.  whereas  the  unaveraged  0.5-in  data  follow 
at  K'st  a  complicated  lognormal  plus  noise  model.  Since 
decimation  leaves  the  underlying  distribution  of  the 
data  unchanged,  it  is  easier  to  provide  a  tractihle  sta¬ 
tistical  model  for  averaged  data  than  for  decimated 
data. 

Why  are  the  high-magnitude  portions  of  <  probability 
densities  lognormal  or  nearly  so'.’  Dissipation  is  log¬ 
normal  or  approximately  so  in  homogeneous  turbu¬ 
lence:  in  the  thermoclinc.  <  can  be  lognormal  within 
individual  turbulent  patches  if  the  sample  length  is  at 
least  three  times  the  Kolmogoroff  length  but  much 
smaller  than  the  height  of  the  patch  ('t'ama/aki  and 
Lueck  1990).  Most  ensembles  of  thermoclinc  data, 
however,  consist  of  records  taken  minutes  to  days  apart 
over  distances  of  tens  to  hundreds  of  meters  in  the 
vertical  and  hundreds  of  meters  to  tens  of  kilometers 
in  the  hori/ontal.  Consequently,  the  samples  come 
from  places  with  no  turbulence  as  well  as  from  many 
different  turbulent  patches.  Furthermore,  the  patches 
are  probably  produced  by  a  variety  of  mechanisms  and 
are  observed  at  varying  stages  in  their  life  cycles.  As  a 
result,  these  bulk  ensembles  do  not  satisfy  the  a.ssump- 
tions  made  in  deriving  lognormal  distnbutions  for  ho¬ 
mogeneous  turbulence  (Gurvich  and  ^  aglom  1967) 
or  within  individual  patches  in  the  thermoclinc  { A'a- 
mazaki  and  l.ucck  1990).  As  also  noted  for  X  (Gregg 
19X0 )  as  well  as  for  <  ( Gibson  1981 ).  the  hulk  statistics 
of  X  and  (  are  likely  to  be  approximately  lognormal 
simply  because  they  have  multiplicative  rather  than 
additive  probabilities.  Based  on  the  similarity  of  the 
high-ntagnitude  portions  of  probability  densities  of  < 
and  .Vto.  we  suggest  that  lognormahty  of  <  may  be  an 
approximate  condition  resulting  from  the  production 
of  turbulence  in  proportion  to  skewed  higher  moments 
of  internal  wave  shear.  A  rigorous  lognormal  model 
may  be  possible,  but  wc  do  not  know  enough  to  con¬ 
struct  one. 

Because  most  turbulence  in  the  thermoclinc  is  pro¬ 
duced  by  breaking  internal  waves,  the  statistics  of  bulk 
ensembles  of  (  must  reflect  the  character  of  internal 
waves  as  well  as  the  inherent  intermittence  of  turbu¬ 
lence.  Whether  or  not  the  fourth  shear  moment  is  the 
best  nuHici,  i  depends  on  a  skewed  shear  momeiit  that 
resembles  some  characteristics  of  the  lognormal  dis¬ 
tribution.  With  only  two  datasets,  we  obviously  have 
not  proven  this  relationship,  but  consider  it  a  hypoth¬ 
esis  to  be  tested. 

If  these  results  arc  substantiated  by  further  obser¬ 
vations,  they  have  several  consequences,  first,  the 
spread  of  the  distribution  of  .S  ;,,  ( or  some  other  higher 
moment )  should  be  an  approximate  upper  bound  for 
die  spread  ol  the  t  distribution.  It  is  an  upiXT  bound 
because  the  low-magnitude  pvirtion  of  the  shear  dis¬ 
tribution  on  average  vlocs  not  produce  turbulence  or 
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it  generates  only  low  t's.  Consistent  with  this,  for  the 
independent  lO-m  samples  we  observe  ff|,„  =  1. 1 9  for 
PATCHEX  and  1.49  for  PATCHEX  north,  signifi¬ 
cantly  less  than  =  2.57.  The  spread  for 

P.ATCHEX  north  is  the  larger  because  the  dissipations 
are  larger  and  less  affected  by  noise.  Comparing  shear 
measurements  at  12  sites  in  the  Atlantic  and  Pacific. 
Smart  ( 19K8)  reports  ^  1-10  times  GV176.  As  Atu 
X  /:,  the  more  energetic  of  these  sites  .should  test  our 
hypothesis  that  2.57  is  an  upper  bound  for  fji„;  pro¬ 
duced  by  internal  wave  shear.  Second,  if  lognormality 
is  only  an  approximate  condition,  lognormal  statistics 
should  not  be  relied  on  a  priori,  as  datasets  may  ap¬ 
proach  lognormality  to  varying  degrees.  Third,  the  sta¬ 
tistics  of  dissipation  arc  likely  to  differ  considerably 
when  turbulent  production  results  from  processes  other 
than  random  internal  wave  fields.  For  example,  our 
development  does  not  apply  to  turbulence  produced 
by  hydraulic  jumps  or  solitary  waves,  such  as  those  in 
the  Strait  of  Gibraltar. 

As  a  consequence  of  the  dependence  of  microstruc¬ 
ture  statistics  on  the  internal  wave  field,  microstructure 
measurements  from  different  times  and  places  cannot 
be  combined  into  a  single  ensemble  for  statistical  anal¬ 
ysis.  For  example.  Gregg  ( 1977)  compares  three  sets 
of  Cox  numbers  obtained  during  different  seasons  at 
28'’N.  !55°W  and  notes  that  their  variability  suggests 
modulation  of  microstructure  intensity  by  mesoscale 
eddies  or  seasonal  effects.  Subsequently,  Briscoe  and 
Weller  ( 1984)  reported  an  apparent  yearly  cycle  in  the 
energy  of  high-frequency  internal  waves  in  which  /; 
varies  by  factors  of  2-1  about  the  long-term  mean. 
Nevertheless,  Baker  and  Gibson  (  1987)  combined  the 
three  sets  of  Cox  numbers  into  a  single  ensemble  with¬ 
out  considering  whether  the  three  sets  were  indepen¬ 
dent.  Proper  treatment  would  consider  each  cruise  as 
an  ensemble  with  its  own  statistics  and  then  examine 
the  statistics  of  the  grand  average. 

How  much  sampling  is  optimum  for  The  27 
PATCHEX  profiles  give  T:,,.  =  0.9.17  and  =  107 
compared  with  =  0.82.1  and  t  -  1-25  for  the 
5  PATCHEX  north  profiles.  Because  we  do  not  believe 
our  calibrations  are  better  than  the  P.A  fCHEX  north 
confidence  limits,  the  additional  sampling  during 
PATCHEX  did  not  improve  our  overall  confidence  in 
the  averages.  If  our  objective  had  been  to  estimate 
whether  a  large  area  had  a  vertical  edtiy  diffusivity  of 
10  "or  10  '  m’ s  '.  a  more  accurate  result  wt)uld  have 
been  obtained  by  sampling  to  obtain  *  lOO'f  accuracy 
at  each  of  many  sites.  The  sites  should  be  separated  by 
about  20  km.  the  approximate  correlation  distance  for 
internal  wave  statistics  (D’Asaro  and  Perkins  1984). 
As.suming  fT|n,  1 .2-2.57  requires .11-282.  equiv¬ 
alent  to  1-4  MSP  profiles  per  site  for  ±  I00'7  accuracy, 
If(T|„;  1.5.  specifying  *  10';  accuracy  for  95'/  con¬ 

fidence  limits,  as  in  the  abstract  of  Baker  and  Gibson 
(  1987  ),  increases  sampling  requirements  from  n  60 


to  //  2000.  or  25  MSP  profiles  pei  site  mstrad  of 

only  one. 

.1ikniwh-Ji’i>h‘nis.  The  Office  ol  Naval  Research 
funded  collection  of  these  data  under  ( Ontraci 
N00014-84-C’-01  I  I  and  analysis  under  Contract 
N00014-86-K-0690.  which  suppvrrtcd  the  .Mixing  to 
Mesoscale  University  Research  Initiative.  11.  ■)  ama/aki 
and  D.  Winkel  gave  valuable  comments,  and  we  arc 
indebted  to  the  reviewer  of  an  earlier  version  of  the 
paper  whose  comment  "Statistical  bafflegab  is  no  sub¬ 
stitute  for  physics  .  .  .  garbage  in,  garbage  out  as  they 
say”  strengthened  our  resolve  to  improv  e  the  statistical 
training  of  oceanographers. 

Al’t'l  \l)l\  A 

Details  on  Statistical  l  ests 

a.  I  erfical  uncorrelaiedncss  of  shear  eoiniymenls 

To  examine  the  vertical  correlation  structure  of  .S\ 
or  .V, .  we  apply  the  run  lest  and  the  spectral  distribution 
function  separately.  The  run  test  is  a  simple  nonpara- 
metric  test  of  uncorrelatedness  ( Bendat  and  Piersol 
1971  ).  A  run  is  a  sequence  of  consecutive  values  lying 
on  one  side  of  the  median.  If  the  n  values  in  the  scries 
are  random  and  independent,  the  number  of  runs.  R. 
is  a  random  variable  with  an  expected  value 

/•;(/?)  =  (2i;./t  //;) -k  1  (Al) 

and  a  variance 

(xj;  =  [2/;./t  {2n.n  -  n)]/[/r(n  -  I)].  (.A2) 

where  n,  and  n  arc  the  number  of  values  greater  and 
less  than  the  median,  and  -i  /i  .  For  large 

/•-'(/?)  =  n/2  and  (Jk  ^  ^n/4.  The  hypothesis  of  un- 
corrclatcdness  can  be  rejected  at  the  5'7  level  of  sig¬ 
nificance  if  R  is  outside  the  interval 

(/-.■(/?)  hiR)  \  I.96.T/,).  (Al) 

For  instance,  each  profile  contains  81. S',  samples  be¬ 
tween  I  and  9  MPa.  leading  us  to  expect  40.5  runs  in 
each  record  with  95'')  confidence  limits  off  12.2,  49.8 ), 

We  applied  the  run  test  to  ail  .V,  profiles  and  to  all 
.V,  profiles.  For  PA  ICHEX,  R  lies  outside  the  95' < 
confidence  limits  in  three  .S\  and  four  .S',  profiles  (  I  able 
A I  ).  For  P.A'K'Hli.X  north,  on  the  other  hand,  one  S', 
and  three  .S',  profiles  are  out  of  bounds.  F'or  comparison, 
with  a  level  of  significance  of  5'  ;  and  with  truly  un¬ 
correlated  data,  we  expect  incorrect  rejections  in  ap¬ 
proximately  0.05  >  27  1,15  profiles  out  of  27,  The 

number  of  incorrect  rejections  itself  varies  .ibcuit  Ms 
expected  value,  fo  assess  the  rnagnitiKlc  of  this  vari¬ 
ation,  we  simulatetl  normally  distributed  white  noise 
with  zero  mean  and  the  same  standard  deviation  as 
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Table  Al.  Tests  for  independence.  The  run  test  (RT)  fora  profile 
is  rejected  as  containing  independent  samples  when  R,  the  number 
of  runs,  exceeds  the  95%  confidence  limits  (A3).  The  sptectral  distri¬ 
bution  function  (SDF)  of  a  profile  is  rejected  when  any  point  exceeds 
the  95%  confidence  limits  for  the  Kolr.iogoroff-Smirnotf  test.  The 
number  of  rejections  is  left  of  the  solidus,  and  the  number  of  runs  is 
right.  Only  for  S,  and  S,  for  PATCHE.X  north  is  (he  same  profile 
(MSP  120)  rejected  by  both  tests.  Each  profile  contains  XI  samples 
of  1 0-m  shear. 


Rejected 
by  RT 

Rejected 
by  SDF 

Rejected 
by  both 

PATCHEX 

5, 

3/27 

0/27 

0/27 

Sy 

4/27 

2/27 

0/27 

2/n 

2/27 

0/27 

5t„ 

2/27 

1/27 

0/27 

PATCHEX  north 

.5, 

1/5 

2/5 

1/5 

Sy 

3/5 

2/5 

1/5 

0/5 

0/5 

0/5 

S*,o 

0/5 

0/5 

0/5 

the  shear  components.  We  grouped  the  simulated  data 
into  27  profiles  of  81  samples  each  and  applied  the  run 
test  in  the  same  fashion  as  we  did  with  the  real 
PATCHEX  data.  By  repeating  this  procedure  1467 
times,  we  observed  2  or  more  rejections  out  of  27  pro¬ 
files  in  34%  of  the  repetitions;  3  or  more  in  12%  of  the 
repetitions;  and  4  or  more  in  3%  of  the  repetitions  (the 
maximum  number  of  rejections  we  observed  was  6). 
Therefore,  the  number  of  rejections  for  PATCHEX 
does  not  suggest  gross  departures  from  uncorrelatedness 
for  or  for  Sy.  This  is  not  true,  however,  for 
PATCHEX  north;  rejections  exceed  the  number  ex- 
pected  by  chance,  indicating  vertical  correlatedness  in 
both  5,  and  Sy. 

The  second  test  uses  the  spectral  distribution  func¬ 
tion,  which  is  the  cumulative  integral  of  the  spectral 
density  function.  (The  empirical  version  of  this  func¬ 
tion  is  known  as  the  cumulative  periodogram.)  Un¬ 
correlated  samples  have  a  white  ( flat )  spectral  density 
function  and,  equivalently,  a  linear  spectral  distribution 
function.  Consequently,  an  empirical  spectral  distri¬ 
bution  function  of  uncorrelated  data  should  fluctuate 
about  a  straight  line,  and  this  is  used  to  test  for  uncor¬ 
relatedness  (Jenkins  and  Watts  1969).  A  profile  is  re¬ 
jected  as  being  white  if  any  point  of  its  spectral  distri¬ 
bution  function  departs  from  the  expected  straight  line 
by  more  than  the  95%  confidence  limits  using  the  Kol- 
mogoroff-Smimoff  test  ( Fig.  A 1 ).  As  is  evident  in  Ta¬ 
ble  A I ,  the  number  of  rejections  is  similar  to  that  for 
the  run  test,  but  the  pattern  of  rejections  is  random. 
Only  one  Sx  and  one  Sy  profile  are  rejected  by  both 
tests,  which  we  take  as  further  evidence  for  no  gross 
departures  from  uncorrelatedness. 


h.  Uncorrelatedness  between  simultaneous  shear 

components 

Taking  all  pairs  of  and  -S’,  for  each  cruise,  we 
compute  Spearman’s  rank-order  correlation  coefficient 
and  Kendall’s  r.  Both  coefficients  are  well-accepted 
nonparametric  tests  whose  validity  does  not  depend 
on  the  assumption  of  normality  (Conover  1980).  We 
compute  the  coefficients  using  numerical  routines  of 
Press  et  al.  ( 1988 ).  In  all  cases  we  cannot  reject  ( at  the 
0.05  level  of  significance)  the  null  hypothesis  that  .S’, 
and  Sy  are  uncorrelated.  Neither  test  comes  close  to 
rejecting  the  null  hypothesis;  the  closest  is  Kendall’s  t 
=  0.0133  for  PATCHEX.  which  can  be  rejected  only 
at  the  35%)  level  of  significance.  For  comparison,  t 
=  —0.0067  for  PATCHEX  north  can  be  rejected  only 
with  an  84%;  level  of  significance.  Therefore,  we  can 
treat  5,  and  S,  as  uncorrelated  with  each  other. 

c.  Lilliefors  statistics  for  95%  confidence  limits 

on  q-q  plots 

The  95%  confidence  limits  are  derived  from  Lilliefors 
statistics  for  normal  or  exponential  distributions  ( Ma¬ 
son  and  Bell  1986).  The  Lilliefors  statistic  is  a  modi¬ 
fication  of  the  Kolmogoroff-Smimoff  st  distic  to  com¬ 
pensate  for  using  estimated  parameters.  Computed  by 
Monte  Carlo  simulations,  the  modification  generally 
gives  a  smaller  acceptance  region  when  parameters 
must  be  estimated  than  when  they  are  assumed  to  be 
known.  Mason  and  Bell  extend  previous  calculations 
by  Stephens  ( 1974)  for  estimated  parameters  and  revise 
Stephens’  critica[  values  using  improved  Monte  Carlo 
techniques.  For  Sx  and  5,.  in  Fig.  2,  we  assume  a  known 
mean  of  zero,  corresponding  to  normal  case  2  of  Mason 
and  Bell.  For  q-q  plots  of  Sio  and  we  use  the  ex- 


Ficj.  .M.  Spectral  distribution  functions  (aka  cumulative  penod- 
ograms)  for  PATCHEX  north  .5,  profiles.  Cumulative  variances, 
normaliz  xl  to  maxima  of  one.  arc  plotted  on  the  y  axis  versus  vertical 
wavenun  her  on  the  .v  axis.  Spectral  distribution  functions  of  white 
noise  would  plot  along  the  diagonal.  The  dashed  lines  arc  95%  con¬ 
fidence  limits  of  the  modified  KolmogorofT-Smirnofl'  test. 
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ponential  case.  For  t,  we  must  estimate  both  the  mean 
and  standard  deviation  from  the  observed  distribution 
and  use  case  3  of  Mason  and  Bell. 


APPENDIX  B 

Derivation  of  Probability  Densities 

a.  Probability  density  of  Sw 

Rather  than  comparing  the  probability  density  of 
5?o  with  observahons,  it  is  more  convenient  to  use  z 
*  logi{)5?o  =  ln5?o/lnlO.  Following  standard  prcxe- 
dures  (Papoulis  1984),  we  take  y  ^  ^10  and 


G(y)  =  z  =  logioy  =  logio-Sfo  =  In^io/lnlO. 

It  has  one  root,  yj  =  e"' and  the  derivative  of  the 
function  is  G'{y)  =  1  /{y  InlO).  Consequently, 

„  _  P,(y,)  In lOf 

b.  Probability  density  of  Sw 

Following  the  same  procedures,  we  return  to  (8)  to 
obtain  the  probability  density  of  .Sto.  Defining  w 
—  ^(>')  =  y  has  only  one  positive  root,  yi  =  Vw. 
Therefore, 

fory*0. 

k(yi)l  4y(rJ 

H'-y^  =  Sto-  (B2) 

The  first  and  second  moments  are 


wP^dw  =  ScTj 


E[w'^] 


f  w^P^dw  = 

J  •'  cc 


384fff.  (B4) 


Therefore,  the  variance  of  is  also  a  simple  function 
of  <Tj, 

al  =  E[n'"]  -  (EIvvI)=  =  320(7?,  (B5) 

and  the  standard  deviation  is  always  2.2  times  the 
mean,  =  17.9<t?  =  2.2£[vv]. 

For  later  use  with  lognormal  statistics,  we  define  z 
^  h{w)  =  Inw,  which  has  one  root,  m’i  =  c*.  Then 


and  the  first  moment  is 


E[z]  =  ln(4a?)  -  2y, 


where  7  =  0.5772 1  •  •  •  is  the  Euler  constant.  The  sec¬ 
ond  moment  is 

£[zM  =  ln'(4(7?)  -  47  ln(4(7?)  +  4^^  +  7'! 


yielding 


ol  =  (2/3)T^ 


Oi  =  <r,aSto  =  2.57. 


(BIO) 


Therefore,  is  constant,  independent  of  variations 
in  ffj. 

To  compare  this  with  observed  distributions,  we  plot 

= . 

4(7? 

zz  =  logioW  =  logio5to.  (Bll) 

Using  (4)  to  scale  S*q  as  hw,  the  logarithmic  form 


P,  = - -r-zi - .  t“logio«m.  (B12) 

4ccr: 


APPENDIX  C 

Statistics  of  .Sio  and  ,$10 

a.  Vertical  uncorrelaiedness  of  Sw  ond  Sw 

No  PATCHEX  north  profiles  are  rejected  by  the 
run  test  or  by  the  spectral  distribution  function  (Table 
Al ),  and  only  2  of  the  27  PATCHEX  profiles  are  re¬ 
jected,  none  by  both  tests.  Consequently,  these  tests 
give  no  indication  of  significant  correlatedness  in  Efo 
and  Eto-  Presumably,  the  moments  have  fewer  rejec¬ 
tions  than  the  shear  components  because  they  depend 
on  both  Sx  and  5,.,  which  are  uncorrelated  with  each 
other. 


b.  Temporal  correlatedness 

To  examine  temporal  correlation  for  PATCHEX. 
we  average  5,  ,  5,  ,  Sfo,  and  5?o  over  100-m  segments 
chosen  to  minimize  decorrelation  due  to  vertical  dis¬ 
placements  by  internal  waves  and  the  internal  tide. 
This  results  in  eight  time  series  for  each  parameter, 
with  each  set  containing  27  values,  one  for  each  profile. 
None  of  the  run  tests  for  £,  or  5,.  must  be  rejected, 
but  two  of  the  eight  series  for  5,0  and  Sw  must  be 
rejected. 

c.  Decrease  in  variances  with  averaging 

The  probability  densities  of  vS|o  and  Sw  allow  us  to 
predict  how  arithmetic  averaging  reduces  the  variances 
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1 ABLK  Cl.  IX'crease  in  o  with  a\craging  length,  /  In  each  enlumn  the  ohsenetl  wiliiev  ol  .'.,,4;,,  01  are  gi\en  lust  .nut  Uilliiwei,l  in 
parentheses  by  the  analytic  \alues  determined  Irom  (Ci!).  f(  3).  and  ((  4l.  .Also  ssithin  the  paienllKses  are  standard  des  lalmiis  oi  die  esliinale'. 
as  determined  I'rom  Monte  Carlo  simulations. 


P.MCHKX 

P\U  111:' 

X  r.unh 

(m) 

10 

1.27 

(1.28  ± 

0.0.3) 

2.55 

(2.5h  ± 

0.05) 

1.18  (1,28 

•  0.07) 

2,  '(>  (2  5(1 

•  0  1  3| 

UK) 

0.42 

(0..32  ± 

0.02) 

0.85 

(0.68  A 

0.0.3) 

0..32  ((I..12  : 

■  0.04) 

0.(i2  (()  (>'' 

'  (JOKi 

8(K) 

0.19 

(0.11  ± 

0.02) 

0..3.3 

(0.24  + 

0.05) 

0.17  (0.1  1  ; 

'  ().0.3| 

0.,)3  (0  21 

■  0081 

of  these  moments.  Owing  to  the  large  skewness  of  the 
probability  densities,  we  describe  the  variances  with 
ffi’nsio  and  which  also  facilitates  comparing  shear 
variances  with  variances  of  t. 

If  the  shear  components  are  independent,  5h)  is  dis¬ 
tributed  as  an  X‘  random  variable  with  two  degrees  of 
freedom;  that  is, 

•^Tn.i  —  cXi,  (Cl  ) 

where  ^lo,,  represents  the  10-m  shear  variances  for  / 
=  1-2187,  and  5h>,,  and  Sj,)./  are  independent  for  1 
#  7.  Then 

=  var(lnST„,,)  =  var(ln(<X5))  =  var(tnX5). 

(C2) 


That  is,  the  proportionality  factor,  c,  does  not  aft'ect 
the  variance. 

These  variances  are  given  analytically  ( Bartlett  and 
Kendal!  1946)  as 

vardnX?/^)  =  i7'(  A),  (C3) 


where  ' '( A')  is  the  derivative  of  the  digamma  function 
(sometimes  called  the  Irigamma  function  ).  F'or  A'  =  I 
we  have  i^'(  I  )  -  7r'/6,  yielding  ffin.v;,,  =  tt/ Vb  =  1.28. 
This  is  in  good  agreement  with  the  observed 
PATCHEX  data  and  within  10'7  of  the  observed 
PATCHF.X  north  data  (Table  Cl).  We  use  Monte 
Carlo  methods  to  assess  the  expected  variability  in 
estimates  of  var  [  inx  ] . 

When  averaging  .S  |,i  ,  over  lOO  and  800  m,  we  arc 
adding  independent  scaled  X“  random  variables  with 
two  degrees  of  freedom,  yielding  scaled  X"  random 
variables  with  more  degrees  of  freedom.  Therefore,  we 
have 


var(  InX  50) 


vardnxrwi).  (C4) 


For  the  lOO-m  PATCHEX  averages  and  for  all  of 
the  800-m  averages,  the  observed  variances  exceed  the 
predictions  (Table  Cl  ) — in  some  cases  as  much  as 


50'-;-.  The  discrepancies  arc  largest  for  the  800-m  av¬ 
erages.  suggesting  that  a  weak  temporal  coherence  may 
be  partly  responsible.  In  view  of  the  discrepancies,  the 
increase  in  degrees  of  freedom  with  averaging  is  not  as 
large  as  expected. 


APPKNDI.X  D 

.Vssessing  the  Lognormal  Plus  Noise  Model 
We  investigate  the  hypothesis  that  the  ().5-m  uncor- 
rccted  t  data  in  a  particular  depth  range  can  be  regarded 
as  the  summation  of(  1  )  a  lognormally  distributed  ran¬ 
dom  variable  with  parameters  and  (7|,„  and  ( 2  )  ran¬ 
dom  noise.  Specifically,  we  assume  that  e  -L  .V. 
where  A  is  a  lognormal  random  variable  with 

Al  Aj  =  c" and  var  '  A  J 

-  A{  A-}  -  (A|  Al)- -  <'-"'‘'’^(0'’^  I).  (Di) 

That  is.  In(.V)  is  Gau.ssian  distributed  with  mean 
and  variance  .  The  noise  variate  .V  is  assumed  to  be 
independent  of  .V  and  to  follow  the  observed  distri¬ 
bution  from  TROPIC'  HliAT  2  between  9  and  9.5  MPa. 
By  ordering  .1/  ^  1480  noise  values  as  «,  <  «-  <  •  •  • 
<  «>/and  assuming  each  equally  likely,  the  probability 
of  a  value  less  than  ,  is 

PI'V  -  '6]  -  (02) 

We  next  determine  the  cumulative  distribution 
function  for  observations  of  dissipation  plus  noise. 

/■A  \  )  P|«  v]  -  P[A  4  .Y  a] 

M 

“  ^  ■  A  and  .\  fh] 

!  I 
M 

P[  A  +  A  ^  vl  .v  /!,]  P[  .v  n,] 

I  1 

1 

-  P|  A  ^  V  n,] 

1  '' 

1  /'d  '  'O-  (0.0 
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where  Fa(  • )  is  the  cumulative  distribution  function  of 
the  lognormally  distributed  variate  A;  that  is. 


F\{  X  -  nj)  = 


ln(x  -  rij)  -  ma 


if  .V  >  n, 

otherwise, 

(D4) 


and  4>(  • )  is  the  cumulative  distribution  function  for  a 
Gaussian  random  variable  with  zero  mean  and  unit 
variance. 

Given  a  sample  «i ,  <2,  of  <  measurements  in 

a  particular  depth  range,  we  fit  our  model  using  a  min¬ 
imum  chi-square  procedure  (Conover  1980,  pp.  195- 
197).  We  use  K  bins  bounded  by  x, .  1  and  x,  such  that 

0  =  A'o  <  X|  <  .X’2  <  •  •  •  <  X/i- I  <  Xa  =  00  .  ( D5) 

The  /cth  class  of  observations  for  /:  =  1  up  to  K  is 
defined  to  be  the  set  of  all  t  measurements  satisfying 
Xa  -i  6/  <  xa.  Let  Ok  be  the  number  of  <  measurements 
in  the  Arth  class  and  let  Ek{^l^,  <Ta)  be  the  expected 
number  of  measurements  in  that  class  under  the  as¬ 
sumption  that  our  proposed  model  is  true  with  log¬ 
normal  parameters  and  <t  v .  We  then  have 


Ekit^x.  (Tx)  -  Lx  (E,{Xk)  -  E,{Xk~i)),  (06) 


which  can  be  computed  for  a  given  and  (Ta  as  out¬ 
lined.  We  can  then  fonn  the  statistic 


T(nx,  a-a)  = 


2  (Ok  -  Lk(.Hx^<r,)f 
A.;  Ek(^lx<<rx) 


2  — — 
A-/  Ek(fkx<  <rx) 


-  L, 


(D7) 


where  we  have  made  use  of  the  fact  that 


/  dEx 

.(  -Va 

-  rij)  dEx\ 

1  Va  I 

\  ' 

dti 

\ 

dfi. 

IfO  <  Xa  - 

-  n,  < 

X  , 

wc  have 

(7<h( 

'ln(.VA  -  «;) 

-  Ua 

a(  Va  - 

'h)  .. 

\ 

,  <^A 

1  /ln{  VA  -  -  Ma 

- ^  — - - 

CTa  \  'fx 

where  </)(x)  =  c  '"'‘/l/23r  (DIO) 

On  the  other  hand,  if  either  va  ^  h,  or  Xk  ~  a  ,  wc  have 
^Ca(-Va  -  n,)/dft  ^  --  0.  Likewise,  we  have 

Ox)  _  I  y  Ol  _ 

dffA  .\f  1^1  Eiikix- ‘^x) 


M 


;  I 


/0/-\(  a'a  -  n,) 

V  <^<7  A 


dExiXk  I  -  n,) 
d(Tx 


(D1  1  ) 


If  0  <  Xa  -  n,  <  X  .  we  have 
dEx(  Xk  -  n,) 
dffx 

ln(.VA  -  n,)  J\n{Xk  -  n,)  -  ma 

= - — -  0 - 

(Ta  \  Ox 

otherwise,  we  have  dE^iXk  -  n,)ldax  =  0. 

Typical  fits  are  shown  in  Fig.  Dl.  and  the  resulting 
estimates  for  ua  and  are  shown  in  Fig.  7.  The 
corresponding  values  for  7(^a,  hx)  are  given  in  Ta¬ 
ble  Dl. 

If  the  </'s  were  independent  of  each  other,  w  e  could 
assess  whether  the  observed  departures  from  the  fitted 


(D12) 


K  K 

2  Ga  =  Z  Ekinx.  Ox)  = 


A- 1  kl 


L. 


(D8) 


The  minimum  chi-square  procedure  is  to  estimate  nx 
and  Ox  as  those  values  such  that  /'(ua-  'Ta)  is  mini¬ 
mized. 

To  find  Aa  and  ox  ( here  the  tilde  indicates  a  sample 
estimate),  we  must  use  a  nonlinear  optimization  rou¬ 
tine,  for  which  it  is  helpful  to  have  expressions  for  the 
partial  derivatives  of  1  inx-  o x)  with  respect  to  a/\  and 
fTA-  Merc 

OTjflx,  Ox)  ^  y  Oj _  ^Ek(klx~  Ox) 

^fix  k  i^'^kifkX^Ox)  (hx 

^  o'  l^f-'-(Xk)  r9/-;(XA  |)\ 

^  lEji^s^Ox)\  OfXx  (ifix  j 


M 


z 

A-  I 


Oi 


Eli  tlx.  Ox) 


M 

z 

I  I 


t  Ui.  1)1.  l  Aanipk's  (if  v’  fils  to  hislograms  ol  ohsorscU 
l).5-m  log, 1,1  taken  m  bins  of  I  MPa. 
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Table  D1.  Chi-square  statistics  forfitsto0.5-m  «:  L  is  the  number  of  samples  in  each  1-MPa  interval  used  for  fitting,  and  K  is  the  number 
of  classes  into  which  the  observations  were  grouped  (this  number  varies  between  28  and  38  to  ensure  that  each  class  contains  at  least  10 
observations).  If  the  L  samples  were  independently  distributed,  the  chi-square  statistic  would  approximately  follow  a  chi-square 

distribution  with  -  3  degrees  of  freedom  (here  and  i*  are  the  minimum  chi-square  estimates  of  the  lognormal  parameters  ^  and  a  ,i. 
Because  the  samples  are  correlated,  assessment  of  is  problematic,  but  a  Monte  Carlo  study  of  normally  distributed  band-limited 

white  noise  indicates  that  we  have  reasonable  goodness  of  fit  for  PATCHEX  north,  whereas  the  opposite  is  true  for  PATCHEX 


MPa 

PATCHEX 

PATCHEX  north 

L 

K 

71ma,«a) 

L 

K 

1-2 

5361 

38 

332.1 

993 

36 

78.4 

2-3 

5372 

35 

249.2 

996 

33 

67.3 

3-4 

5382 

32 

326.8 

997 

33 

107,9 

4-5 

5382 

34 

305.9 

1000 

37 

53.0 

5-6 

5384 

32 

299.2 

996 

30 

51.4 

6-7 

5396 

29 

292.5 

998 

38 

80.8 

7-8 

5389 

28 

398.0 

998 

33 

52.9 

8-9 

5393 

30 

316.1 

997 

36 

48.9 

distributions 

are  consistent 

with  random 

error  by  Let  jCp 

be  the  pth  quantile  of  the 

distribution. 

comparing  the  value  of  Tips 

,  5a)  to  a  chi-square  ran-  Then 

dom  variable  with  K—'i  degrees  of  freedom  (Conover 
1980).  Unfortunately,  the  «/  are  not  anywhere  close  to 
being  independent,  but  based  on  the  assumption  that 
correlation  and  distributional  properties  can  be  decou' 
pled  in  a  chi-square  statistic,  we  can  crudely  determine 
the  effect  of  different  degrees  of  dependence  on  the  chi- 
square  statistic  by  the  following  procedure. 

For  an  integer  A/  >  1 ,  consider  a  normally  distrib¬ 
uted  band-limited  while  noise  process  with  a  low-fre¬ 
quency  cutoff  of  1  /  2  M .  The  theoretical  autocovariance 
sequence  for  such  noise  is  positive  for  lags  less  than  M 
and  is  zero  at  lag  M\  hence,  by  choosing  M  in  the  range 
from  10  to  20,  the  autocovariance  sequence  for  the 
band-limited  process  will  mimic  the  observed  sequence 
for  the  0.5-m  t/,  which  damps  down  to  zero  typically 
in  the  range  from  5  to  10  m.  We  then  generate  samples 
from  this  process  with  a  number  of  points  comparable 
to  those  obtained  in  either  PATCHEX  (5380)  or 
PATCHEX  north  ( 1000)  and  fit  the  data  to  a  normal 
distribution  using  the  minimum  chi-square  procedure 
outlined  above.  By  replicating  this  procedure  many  dif¬ 
ferent  times  and  looking  at  the  observed  values  of  the 
chi-square  statistic,  we  find  that  the  range  of  observed 
chi-square  statistics  in  the  simulations  matches  very 
closely  the  range  of  values  for  T(Jix,  as)  shown  in  Table 
D1  for  PATCHEX  north.  For  PATCHEX,  however, 
the  values  are  about  a  factor  of  2  to  3  larger  than  those 
in  the  simulation  study.  Thus,  we  conclude  that  our 
operational  model  fits  the  PATCHEX  north  data  rea¬ 
sonably  well,  while  there  is  evidence  of  a  lack  of  fit  in 
PATCHEX  (this  can  be  seen  in  Fig.  Dl,  where  the 
fitted  distribution  and  the  histogram  disagree  markedly 
in  the  central  portion). 

APPENDIX  E 

Theoretical  q-q  Plots 


P  =  PIc5to.,  <  Xp]  =  P[C<TtX2  Xp] 

=  P[xl^)/l^-/aU  =  P[xl^x'p].  (El) 

where  x'p  is  the  pth  quantile  of  the  X  2  distribution.  We 
thus  have 

Xp=(C(r^Xp)^  (E2) 

Next,  let  Tp  be  the  pth  quantile  of  the  e  disiril/ulion. 


g-q  plots  of  lO-m  t  and 


FiCi.  El.  Plots  of  theoretical  Q-q  distributions  for  logio«„  and 
log,o<  agree  well  at  high  magnitudes  hut  not  at  low  magnitud'*s 
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and  let  Z  be  a  noniial  random  variable  with  zero  mean 
and  unit  variance;  then 

log(a)  -  Mm;  ^  logl^p)  -  Min;' 

^Iii<  <^lnl 


p  =  p[a  <  vp]  =  p 


=  p 


10g(>V)  -  Min; 


^ln« 


-P[Z^y',], 


(E3) 


where  y'p  is  the  pth  quantile  of  the  standard  normal 
distribution.  We  thus  have 


>’p  =  exp(  Min;  +  y'p  a,„:) .  ( E4 ) 

The  q~q  plots  of  these  two  theoretical  distributions 
agree  well  at  high  magnitudes  (Fig.  El ). 

APPENDIX  F 

Coofidcacc  Limits  for  and  £(5io) 


[LnX<5fo>,  C''nX<5?0>],  (F7) 

where 

In  In 

-  and  - .  (F8) 

^2b.(I+p>/2  X2n,(l-p)/2 

To  evaluate  X\„^p,  we  use  Eq.  (3)  of  Zar  ( 1978);  to 
input  X,  the  standardized  normal  deviate,  we  use  the 
approximation  for  4>”'(p)  in  Table  6.5  of  Chambers 
etal.  (1983). 

b.  Confidence  limits  for  E(S%) 

As  noted  in  the  main  text,  confidence  limits  for 
E(S*o)  can  be  based  either  on  (5io)  or  on 

(Sto)  ^  -  Z  5to.,.  (F9) 

"  ,-1 


a.  Confidence  limits  for  E{S]o) 

Previously  we  found  that  5?o  is  distributed  as  a  re¬ 
scaled  X  ^  random  variable  with  two  degrees  of  freedom: 


S\i 


i  y2 
-  Xj 


(FI) 


using  E(Siq)  =  la].  Averaging  n  of  these  Sio,,  yields 


I- 1 


In 


(F2) 


or,  equivalently. 


2n<£fo>  ^ 

£(5?o) 


xL, 


(F3) 


where  X2„  denotes  a  X^  random  variable  with  In  de¬ 
grees  of  freedom. 

To  obtain  confidence  limits  for  E{S\q),  let  X  2n,p 
represent  the  p  X  100th  percentage  point  of  the  X^ 
distribution  with  In  degrees  of  freedom;  that  is, 

P[xL  ^  Xl„,p]  =  p,  (F4) 


where  P[A]  refers  to  the  probability  of  the  event  A. 
Then 


2  ,  =p.  (F5) 


pm 

lo) 

This  is  equivalent  to 

2n(S]o} 


L^L.(I  +  p)/2 


EiSio) 


2n(S]o} 

X2n,(l-p)/2j 


=  p.  (F6) 


Therefore,  the  p  X  100%  confidence  interval  for 
£(S]o)  is 


Let  Q„  denote  the  average  of  the  squares  of  n  indepen¬ 
dently  distributed  xi  random  variables,  and  let  Q„j, 
represent  the  p  X  100th  percentage  point  of  the  distri¬ 
bution  of  Q„\  that  is,  P[Cn.<  Qnj>]  =  p.  py.  100% 
confidence  interval  for  flSto)  is  given  by 

[L;  X  <5to>.  i/'n  X  <5to>],  (FIO) 

where 

E’n  =  -p-^ -  and  t/;  =  — - .  (Fll) 

tin.(l+P>/2  C^n.()-p)/2 

The  percentage  points  Q„j,  can  be  determined  either 
by  Monte  Carlo  simulations  or,  if  n  3=  100,  by  an  Edge- 
worth  approximation  to  the  distribution  of  Q„  (see  sec¬ 
tion  1 .5.D  of  Bickel  and  Doksum  ( 1977 )] . 


APPENDIX  G 

Comments  on  Baker  and  Gibson  (1987) 

Baker  and  Gibson  ( 1987 )  urge  using  lognormal  sta¬ 
tistics  to  estimate  (e)  in  the  thermocline.  We  agree 
that  lognormal  statistics  are  useful  for  planning  field 
measurements,  but  do  not  concur  that  they  should  be 
routine  for  post  facto  data  analysis.  In  the  first  place, 
not  all  datasets  are  lognormal.  As  noted  in  the  text, 
our  t  profiles  are  definitely  not  lognormal  unless  scaled 
with  <  )  and  corrected  for  noise.  Second,  at  present 

there  is  no  rigorous  model  demonstrating  under  what 
conditions  e  should  be  lognormal  in  the  thermocline. 
Our  hypothesis  in  this  paper  is  that  lognormality  of 
scaled  and  noise-corrected  data  is  an  approximate 
rather  than  an  exact  condition. 

T’urd,  Baker  ar^  Gibson  (1987)  claim  that  the  mle- 
bas'.d  estimator  A'mie  of  the  parameter  p  has  smaller 
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Table  Gl,  Statistics  of  lO-m  dissipation  between  1  and  9  MPa, 
For  small  samples  of  lognormal  distributions,  arithmetic  estimators 
of  the  mean  are  more  efficient  than  maximum  likelihood  estimators; 
that  is.  varTt'a^  ;  <  var!  .\„,|,;  tor  n  •  j  I  .  '  'r  ji  <,  lO  the  efficienes 
is  infinite. 


n 

\  ur  i  -V  mic  1 

var;.\\„; 

Ratio  (eliieieney ) 

to 

41  381  .397  761 

259  1  56 

159  678 

15 

6.31  252 

172  771 

3.65 

20 

153  17.3 

129  578 

I.I8 

21 

128  .3(X) 

1 23  408 

1.04 

22 

109  528 

1  17  798 

0.93 

25 

74  122 

103  662 

0.71 

50 

16  761 

51  831 

0.32 

too 

6  046 

25  916 

0.23 

X 

— 

— 

0.17 

variance  than  the  sample  mean  This  is  true  for 
moderate-to-large  sample  sizes,  but  it  is  nol  true  for 
very  small  samples.  In  fact,  for  samples  sizes  n  such 
that 

n  <  2rT"  +  1  (0 1 ) 

the  variance  of  A'^ie  is  infinite!  This  condition  is  cor¬ 
rectly  stated  by  Baker  and  Gibson  following  their  Eq. 
(Bl  1 ),  but  they  do  not  note  its  practical  implications. 
This  becomes  an  issue  in  their  section  3.  where  they 
represent  an  example  from  Gregg  (1977)  with  <t- 
=  4.41  and  «  =  10.  For  this  value  of  a*,  n  =  10  is  the 
smallest  sample  size  for  which  var{  Anuc}  is  finite.  If 
we  use  their  Eqs.  (BIO)  and  (Bl  1 ).  we  can  compute 
var  { ,Vn,ic } ,  var  { A'am  j  and  their  ratio  for  small  n .  Table 
Gl  shows  that,  for  «  ^  21,  we  should  prefer  A'am  over 
A'mic  as  our  estimator  of  n.  while  the  opposite  holds  for 
n  >  2 1 .  ( Note  that  their  Fig.  B 1  only  shows  the  relative 
efficiency  for  moderate  n,  namely.  50  and  100.) 

If  we  are  dealing  with  correlated  data  rather  than 
uncorrelated  data  as  Baker  and  Gibson  assume,  the 
minimum  sample  size  for  which  .Vmu.  is  preferable  over 
■Vam  will  certainly  be  larger. 

It  should  also  be  noted  that  confidence  intervals 
based  on  mle’s  are  not  necessarily  valid  when  other 
estimates  are  used.  Thus,  Baker  and  Gibson's  claim 
that  mic-based  confidence  intervals  can  be  used  with 
their  graphical  estimates  is  not  correct.  To  be  specific, 
the  claim  that  their  equation  (7)  still  holds  when  I  ) 
the  graphical  cstimati^r  is  used  and  2 )  there  is  substan¬ 
tial  noise  cannot  be  supported  by  appealing  to  the 
noiseless  case  (as  the  authors  do  in  their  diseussion  of 
Fig.  4h ).  As  can  be  seen  from  a  limiting  case  argument 
( i.e.,  all  noise),  there  be  some  dependence  on  the 
number  of  points  included  in  the  graphical  fit.  but  they 
claim  there  is  none.  Furthermore,  their  graphical  es¬ 
timates  of  )i|n;  and  ff|,„  arc  not  maximum  likelihood 
estimates. 
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